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Measurement of at least three independent parameters, for example, mass, radius and spin frequency, of a neutron star 
is probably the only way to understand the nature of its supranuclear core matter. Such a measurement is extremely 
difficult because of various systematic uncertainties. The lack of knowledge of several system parameter values gives 
rise to such systematics. Low mass X-ray binaries, which contain neutron stars, provide a number of methods to 
constrain the stellar parameters. Joint application of these methods has a great potential to significantly reduce the 
j systematic uncertainties, and hence to measure three independent neutron star parameters accurately. Here we review 
the methods based on (1) thermonuclear X-ray bursts; (2) accretion-powered millisecond-period pulsations; (3) kilohertz 
quasi-periodic oscillations; (4) broad relativistic iron lines; (5) quiescent emissions; and (6) binary orbital motions. 
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1. Introduction 

The nature of matter at ultra-high density, i.e., 5 — 10 
times higher than the nuclear density, but at a relatively 
low temperature (e.g., ~ 10 8 K) cannot be probed by 
heavy-nuclei collision experiments or with observations of 
the early universe (Qzel and Psaltis (2009); van Kerk- 
wijk (2004); Blaschkeet al. (2008); Lattimer and Prakash 
(2007) and references therein). This matter exists in the 
neutron star cores, and hence these stars provide natural 
laboratories to study the supranuclear degenerate matter. 
The study of this dense matter is not only important from 
the point of view of the particle physics, but also can be 
useful to understand the core collapse of massive stars, 
the supernova phenomenon, the existence and properties 
of neutron stars, etc. (e.g., van Kerkwijk (2004)). 

But how to observe and understand the deep interior of 
a neutron star? Based on the assumed microscopic prop- 
erties of the core matter, many equation of state (EoS) 
models have been proposed in the literature. These mod- 
els connect the pressure (p) with the total mass-energy 
density (e) in a degenerate condition. Therefore, for a 
nonspinning neutron star and for a given EoS model, the 
stellar global structure can be calculated by solving the 
Tolman-Oppenheimer-Volkoff (TOV) equation: 

dp(r) G [p(r) + e(r)][m(r) + 47rr 3 p(r)/c 2 ] 
dr c 2 r(r — 2Gm(r)/c 2 ) 

where, m(r) — An f Q r dr'r /2 e(r') is the gravitational mass 
(hereafter, mass) inside a sphere of radius r. For a rapidly 
spinning neutron star, Cook et al. (1994) and other au- 
thors have suggested ways to compute the stellar struc- 
ture (see the dashed curves of Figure 1; see also, Datta et 
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al. (1998); Thampan et al. (1999); Bhattacharyya et al. 
(2000, 2001a,b)). The Figure 1 shows that for a given EoS 
model and a known stellar spin frequency the stable struc- 
tures of neutron stars trace a single mass(M )-radius(i?) 
curve. However, such curves for different EoS models can 
intersect each other (Figure 1). Therefore, it is required 
to accurately measure at least three independent global 
parameters of the same neutron star in order to constrain 
the EoS models, and hence to understand the supranu- 
clear core matter. This is also true for the EoS models 
of strange matter (made of u, d and s quarks; Witten 
(1984); Farhi and Jaffe (1984); Bombaci (1997); Bom- 
baci et al. (2000); Bhattacharyya et al. (2001c); Bagchi 
et al. (2006)). However, some neutron star EoS mod- 
els can be ruled out only by mass measurement. Here is 
how. A stable neutron star's mass cannot be greater than 
a maximum value (M maK ; Lattimer and Prakash (2007)) 
for a given EoS model and spin frequency. Therefore, a 
large measured mass value can rule out some "softer" EoS 
models (for example, "NS1" in Figure 1), that have lower 
M max than the "harder" models (for example, "NS4" in 
Figure 1; see also Lattimer and Prakash (2007)). 
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Figure 1: This figure shows the mass-radius (M — R) space of neu- 
tron stars and the curves corresponding to a few representative EoS 
models (§ 1). The solid curves are for nonspinning stars, while the 
dashed curves (corresponding to some of the solid curves) are for a 
stellar spin frequency of 600 Hz. The red curves are for the usual 
neutron star EoS models: NS1 (A18 model of Akmal et al. (1998)); 
NS2 (Baldo et al. , 1997); NS3 (A18+5d+UIX model of Akmal et 
al. (1998)); NS4 (Sahu et al. , 1993); and NS5 (Pandharipande and 
Smith , 1974); while the blue curves are for the strange matter EoS 
models: SSI (bag constant = 90 MeV fm -3 ; Farhi and Jaffe (1984)) 
and SS2 (bag constant = 60 MeV mi -3 ; Farhi and Jaffe (1984); see 
also § 1). The grey region indicates the M — R space of black holes. 
Here we have shown the M > 1 Mq region, because many neutron 
stars in double neutron star binaries have the mass around 1.35 Mq 
(Thorsett and Chakrabarty , 1999), and therefore the accreting neu- 
tron stars in low-mass X-ray binaries should have even larger mass 
values (see § 4). This figure suggests that the measurement of mass, 
radius and spin frequency of the same neutron star is required to 
constrain the EoS models. 



The above discussion give the motivation for an accu- 
rate measurement of the neutron star parameters. In this 
review, we discuss the measurement methods relevant for 
neutron stars in low-mass X-ray binary (LMXB) systems 
(Bhattacharya and van den Heuvel , 1991). An LMXB is 
an old binary stellar system (typical age ~ 10 9 years) with 
a low-mass companion star (< IMq) and a compact star 
rotating around each other. The companion star may be a 
main sequence star, an evolved star or a white dwarf, while 
the compact star may be a neutron star or a black hole. 
In such a system, the companion star fills its Roche lobe, 
i.e., the critical equipotential surface (see Bhattacharya 
and van den Heuvel (1991)), and matter from it flows to- 
wards the compact star (neutron star in our case). This 
matter cannot move radially towards the neutron star, as 
its initial angular momentum is large. Therefore, the ac- 
creted matter slowly gets rid of the angular momentum 
(Shakura and Sunyaev , 1973), follows spiral paths, and 
hence forms an accretion disk. In neutron star LMXBs, 
the accreted matter eventually hits the stellar surface, and 
generates electromagnetic radiation. The accretion disk 
also emits such radiation by viscous dissipation, and both 
these emissions are powered by the gravitational poten- 
tial energy release. The inner part of the accretion disk, 
and the neutron star surface emit X-rays, as their typical 
temperature is ~ 10 7 K. These X-rays are reprocessed by 
the outer part of the accretion disk, which emits in optical 
wavelengths. Moreover, the companion star may be irra- 
diated by these X-rays. The neutron stars in LMXBs are 
expected to be spun up, due to the angular momentum 
transferred via accretion. This is consistent with the ob- 
served stellar spin frequencies (typically ~ 300 — 600 Hz; 
see § 3.1.6; § 3.2.2). The magnetic fields (w 10 7 - 10 9 G) 
of these neutron stars are typically much lower than that 
of the isolated neutron stars and the stars in high-mass X- 
ray binary (HMXB) systems (Bhattacharya and van den 
Heuvel , 1991; Ibragimov and Poutanen , 2009; Psaltis and 
Chakrabarty , 1999). Several groups have worked on this 
problem of stellar magnetic field evolution (Srinivasan et 
al. , 1990; Jahan Miri and Bhattacharya , 1995; Konar et 
al. , 1995; Choudhuri and Konar , 2002; Konar and Choud- 
huri , 2004). Note that the low magnetic field allows the 
accretion disk to extend up to very close to the neutron 
star. This, in combination with the fact that X-rays can 
be detected both from the star and the inner part of the 
disk for many LMXBs, provides an excellent opportunity 
to measure the neutron star parameters. Moreover, the 
signatures of binary orbital motion detected from the op- 
tical observation of the accretion disk and the companion 
star can be useful to measure the neutron star mass. 

LMXBs are observed mostly in the Galactic plane and 
bulge, and in the globular clusters. Some of them have 
also been discovered from the nearby galaxies. According 
to the catalogue of Liu et al. (2007), 187 LMXBs (con- 
taining either a neutron star or a black hole) were known 
from our galaxy, "Large Magellanic Cloud" and "Small 
Magellanic Cloud" up to that time. As indicated from 
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thermonuclear bursts, mass function, etc., most of these 
LMXBs contain a neutron star, rather than a black hole. 
LMXBs emit primarily in X-ray wavelengths, and the op- 
tical to X-ray luminosity ratio is normally less than 0.1 
for them (see Bhattacharya and van den Heuvel (1991)). 
The central X-ray emission from these sources cannot be 
spatially resolved using the current X-ray imaging instru- 
ments, because the distances of even our Galactic LMXBs 
are very large (typically > 1 kpc). Therefore, one needs 
to rely on the spectral and timing properties in order to 
study these sources in X-rays. 

X-rays from space cannot reach the surface of the Earth. 
Therefore, the X-ray detectors and telescopes are sent above 
the atmosphere, typically by satellites. The current pri- 
mary X-ray space missions are (1) NASA's Rossi X-ray 
Timing Explorer (RXTE); (2) ESA's XMM-Newton; (3) 
NASA's Chandra; and (4) JAXA's Suzaku. They pri- 
marily operate in the w 1 - 10 kcV range (i.e., soft X- 
rays), and contain a variety of instruments. The primary 
instrument of RXTE is a set of five proportional coun- 
ters (called "Proportional Counter Array" or PCA) with 
large collecting area and very good time resolution (« 1 
microsecond). Therefore, PCA is an ideal instrument to 
study the fast timing phenomena, such as burst oscillations 
(§ 3.1.6), accretion-powered pulsations (§ 3.2), kilohertz 
quasi-periodic oscillations (§ 3.3), and also the continuum 
energy spectrum (e.g., § 3.1.3). PCA is also the best in- 
strument to study most of the properties of thermonuclear 
bursts (see § 3.1). However, it is not suitable to study the 
spectral lines because of its poor energy resolution. This 
instrument operates in 2 — 60 KeV band, although the 
effective collecting area becomes small above w 20 keV. 
XMM-Newton satellite contains X-ray telescopes, charge- 
coupled devices and high resolution "Reflection Grating 
Spectrometers (RGS)" . Therefore, this satellite, which op- 
erates in w 0.2 — 12 keV range, is ideal for spectral analysis, 
including narrow spectral lines. The Chandra observatory 
contains an X-ray telescope, charge-coupled devices, high 
resolution camera and transmission gratings. Its angu- 
lar resolution (~ subarcsecond) is the best among all the 
past and current X-ray instruments, and it is also ideal 
to detect and study narrow spectral lines (in w 0.1 — 10 
kcV range). XMM-Newton and Chandra are also useful 
to study the quiescent emission of neutron star LMXBs 
(§ 3.5). The Suzaku satellite also has X-ray telescopes and 
charge-coupled devices. Since it contains an additional 
"Hard X-ray Detector" (HXD), it can study the broad- 
band energy spectrum. Therefore, Suzaku is very useful to 
observe the broad relativistic iron emission lines and the 
associated "disk reflection" spectra (see § 3.4). 

Here is the plan of this review. In § 2, we describe 
the aim and the timeliness of this review. In § 3, we 
briefly discuss a few observed properties of neutron star 
LMXBs, how these properties can be used to measure 
neutron star parameters, the relevant systematics and the 
plausible ways to reduce these systematic uncertainties. 
§ 4 mentions why LMXBs are useful to measure the neu- 



tron star parameters, and § 5 includes a summary of this 
review, as well as a brief discussion on future prospects. 

2. Why This Review and Why Now? 

This review aims to describe the major neutron-star- 
parameter-measurement methods for LMXB systems briefly, 
and to be useful as a quick reference for astrophysicists, 
nuclear physicists, particle physicists, and the scientists 
working in fluid mechanics, general theory of relativity, 
and other relevant fields. This is probably the first time 
that all the methods, including a recently established method 
based on broad relativistic spectral iron emission lines (see 
§ 3.4.2), for LMXB systems have exclusively been compiled 
together. Although, this review does not include the tech- 
niques for all types of neutron star systems, the discussed 
methods form a complete set for the purpose of constrain- 
ing EoS models. This is because at least three independent 
parameters of the same neutron star must be measured to 
constrain these models (see § 1). 

We believe that this review is timely, because of sev- 
eral new discoveries and advancement in the field of neu- 
tron star LMXBs made in the past few years. To men- 
tion a few, (1) Chakrabarty et al. (2003) established 
that the burst oscillations do give the neutron star spin 
frequency (see § 3.1.6); (2) the sample size of accretion- 
powered millisecond-period pulsars has increased from 1 
to 12 (§ 3.2.1); (3) Cottam ct al. (2002) reported the only 
plausible observation of the neutron star surface atomic 
spectral lines (see § 3.1.4); (4) Bhattacharyya and Strohmayer 
(2007b) established the inner accretion disk (or relativis- 
tic) origin of the broad iron line from a neutron star LMXB 
for the first time (see § 3.4.1); and (5) Steeghs and Casares 
(2002) detected the companion star of an LMXB based on 
Bowen blend emission lines, and constrained the neutron 
star mass using these lines for the first time (see § 3.6). 

3. Methods 

X-ray emission of an LMXB from a region close to the 
neutron star surface contains information about the stel- 
lar parameters. This information can be extracted from 
the observed spectral and timing properties of these X- 
rays. Here we discuss the properties of the following five 
X-ray phenomena, and the methods for neutron star pa- 
rameter measurement using them: (1) thermonuclear X- 
ray bursts; (2) accretion-powered millisecond-period pul- 
sations; (3) kilohertz quasi-periodic oscillations; (4) broad 
relativistic iron lines; and (5) quiescent emissions. In ad- 
dition, we discuss a mass measurement technique based on 
the orbital motions of the binary components. 

3.1. Thermonuclear X-ray Bursts 

3.1.1. What is a Thermonuclear X-ray Burst? 

Eruptions in X-rays are observed from many neutron 
star LMXB systems every few hours to days (Strohmayer 
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and Bildsten , 2006; Galloway et al. , 2008a). These are 
called type-I X-ray bursts. The observed X-ray intensity 
sharply increases typically by a factor of~10in«0.5 — 5 
seconds, and then decreases relatively slowly in w 10 — 100 
seconds during such a burst (Galloway ct al. (2008a); 
see also Figure 2). The typical energy emitted in a few 
seconds is <~ 10 39 ergs (Strohmayer and Bildsten , 2006). 
These bursts were first discovered in 1970's (Grindlay et 
al. , 1976; Belian et al. , 1976). Soon after this discov- 
ery it was realized that they originate from intermittent 
unstable nuclear burning of accreted matter on the neu- 
tron star surfaces (Joss , 1977; Lamb and Lamb , 1978; 
Strohmayer and Bildsten , 2006). The neutron star sur- 
face origin was supported by the observational fact that 
the burst emission area matched with the expected surface 
area of a neutron star (Swank ct al. , 1977; Hoffman et al. , 
1977a). Note that the emission area can be estimated from 
the fitting of the energy spectra of the bursts. The unsta- 
ble burning is required because, while the nuclear energy 
released per nucleon by fusion is only a few MeV, the grav- 
itational energy released by a nucleon falling on the neu- 
tron star surface from a large distance is GMm nuc \ eon / R 
(M and R are neutron star mass and radius respectively) , 
which is typically <~ 200 McV. Therefore, if the nuclear 
burning were stable, the nuclear energy would be entirely 
lost in the observed gravitational energy, and the burst 
could not be powered by the nuclear energy release. This 
nuclear burning instability is driven by a nuclear energy 
generation rate which is more temperature sensitive than 
the radiative cooling. Such a "thin shell" instability was 
first theoretically discovered by Schwarzschild and Harm 
(1965) for the asymptotic giant branch phase of stellar 
evolution, and was shown to occur in the thin shell of hy- 
drogen and helium on neutron star surfaces (Hansen and 
van Horn , 1975). But what shows that the nuclear burn- 
ing causes the bursts? A primary evidence of it comes 
from the following observational fact. For a data set from 
a given source, the ratio of the total energy content in all 
the bursts to that in the non-burst emission roughly tallies 
with the ratio of the expected nuclear energy release per 
nucleon to the expected gravitational energy release per 
nucleon. This strongly suggests that the accreted matter 
accumulates on the neutron star surface for a long time, 
and then produces a burst via a quick thermonuclear burn- 
ing. Therefore these bursts are known as "thermonuclear 
bursts" (see the reviews Lewin et al. (1993, 1995) and 
Strohmayer and Bildsten (2006) for detailed descriptions 
of these bursts). 

For many neutron star LMXBs, the accreted mate- 
rial consists of mostly hydrogen, some helium and a small 
amount of heavier elements, although for ultra-compact X- 
ray binaries (UCXBs) the donor companion is a hydrogen 
deficient and heavier element rich dwarf star (in't Zand et 
al. (2008) and references therein). The accreted matter, 
which accumulates on the neutron star surface, undergoes 
hydrostatic compression as more and more material piles 
on, and the ignition column density and temperature are 
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Figure 2: This figure shows the X-ray intensity profile of a typical 
thermonuclear X-ray burst observed from an LMXB. The burst light 
curve is for the entire RXTE PCA band, and with 1 sec binsize. 
After the onset of the burst, the X-ray intensity increases sharply by 
more than an order of magnitude compared to the pre-burst intensity 
(powered by the gravitational potential energy), and then decays 
relatively slowly as the neutron star surface cools down (§ 3.1.1). 

reached in a few hours to days (Schatz and Rchm , 2006). 
The burning typically happens at a depth of « 10 m at 
a column density of ~ 10 8 gm cm -2 . The compression 
rate and the ignition depends on the accretion rate, which 
sets four distinct theoretically expected regimes of burning 
(Bildsten , 2000). 

Regime 1: The temperature of the accumulated material 
exceeds 10 7 K in most cases, and hence the hydrogen burns 
via the CNO (carbon-nitrogen-oxygen) cycle instead of the 
pp (proton-proton) cycle. This causes mixed hydrogen 
and helium bursts triggered by the unstable hydrogen igni- 
tion, which happens for the accretion rate per unit neutron 
star surface area rh < 900 gm cm~ 2 sec -1 (Zcno/0-01) 1 / 2 
(Bildsten , 2000). Here Zqno is the mass fraction of CNO. 
Regime 2: the proton capture timescale becomes shorter 
than the subsequent (3 decay lifetimes at higher tempera- 
tures (T > 8 x 10 7 K). Hence, the hydrogen burns via the 
"hot" CNO cycle of Fowler and Hoyle (1965): 

12 C(p, 7 ) 13 N(p, 7 ) 14 0(/?+) 14 N( P , 7 ) 15 0(/3+) 15 N(p, a) 12 C. 

In this case, the hydrogen burning happens in a ther- 
mally stable manner (i.e., without triggering a thermonu- 
clear burst) simultaneously with the accumulation of mat- 
ter. This allows a helium layer to build up below the 
hydrogen layer. For the rh range 900 gm cm -2 sec -1 
(ZcNo/O.Olj 1 / 2 <m< 2xl0 3 gm cm- 2 sec" 1 (Z C no/0.01) 13 / 
the hydrogen entirely burns before this helium is ignited 
(Bildsten , 2000; Strohmayer and Bildsten , 2006). So 
when the helium ignition happens, a short (~ 10 sec) but 
very intense burst occurs by the unstable triple-alpha re- 
action of the pure helium (3a — > 12 C). Such a burst is 
called a "helium burst" . Many of these bursts are so in- 
tense that the local X-ray luminosity in the neutron star 
atmosphere may exceed the Eddington limit (for which the 
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radiative pressure force balances the gravitational force), 
and the photospheric layers may be lifted off the neutron 
star surface. Such bursts are called photospheric radius 
expansion bursts (Strohmayer and Bildsten (2006) and 
references therein). 

Regime 3: hydrogen burns via the "hot" CNO cycle in a 



and references therein) . The typical released energy, recur- 
rence time and decay time of these bursts are ~ 10 42 ergs, 
~ years and ^1 — 3 hours respectively (Cumming and 
Bildsten , 2001; Strohmayer and Brown , 2002; Cooper 
et al. , 2009). Thermally unstable 12 C fusion (Woosley 

and Taam , 1976; Taam and Picklum , 1978; Brown and 
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stable manner for m > 2x 10 gm cm sec" _i (^ C No/0.01) 13/18 Bildstcn , 1998) at a column depth of - 10" gm cm 
(as for the regime 2). But in this case, enough amount of 
unburnt hydrogen remains present at the time of helium ig- 
nition. This is because at this higher accretion rate, the he- 
lium ignition conditions are satisfied much sooner. There- 
fore in this regime, mixed hydrogen and helium bursts are 
triggered by the helium ignition. During such a burst, 
the thermal instability can produce elements beyond the 
iron group (Hanawa et al. , 1983; Hanawa and Fujimoto 
, 1984; Koike et al. , 1999; Schatz et al. , 2001) via the 
rp (rapid-proton) process of Wallace and Woosley (1981). 
This is because the triggering helium flash makes the tem- 
perature high enough to allow the "breakout reaction" 
15 0(a,7) 19 Ne (Wallace and Woosley , 1981; Schatz et al. 
, 1999; Fisker et al. , 2006) from hot CNO cycle to proceed 
faster than the j3 decays. 19 Ne can return to the hot CNO 
cycle via the chain of reactions (Cooper et al. , 2006): 
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Ne(/3+^) 19 F(p, a) 16 0(p, 7 ) 17 F(p, 7 ) 18 Nc(/3 4 >) 18 F(p, a) 15 0. 

But if 19 Ne captures a proton by the reaction 19 Ne(p, 7) 20 Na, 
the nuclear reactions go out of the hot CNO cycle loop, 
and follow the rp process. This process burns hydrogen 
via successive proton capture and (3 decays, and produces 
a large range of heavy nuclei (Schatz et al. , 1999, 2001). 
The long series of (3 decays makes these bursts typically 
much longer (~ 100 sec) than the helium bursts. 
Regime 4- finally at a very high accretion rate (comparable 
to the Eddington limit), the helium burning temperature 
sensitivity becomes weaker than the cooling rate's sensitiv- 
ity (Ayasli and Joss , 1982; Taam et al. , 1996). Therefore 
in this regime, the stable burning sets in, and the ther- 
monuclear bursts do not occur. However, recently Cooper 
and Narayan (2006) have shown that the bursts should 
not occur even if the accretion rate is more than 30% of 
the Eddington accretion rate, in case 15 0(a, 7 ) 19 Ne reac- 
tion rate is lower than usually assumed. This appears to 
agree with observations (Cooper and Narayan , 2006). Fi- 
nally, Cooper and Narayan (2007a) have suggested that 
for rapidly spinning neutron stars, burning instability sets 
in (i.e., bursts ignite) preferentially near the equator at low 
values of accretion rate and off of the equator at higher ac- 
cretion rates. 

Apart from usual type-I X-ray bursts, bursts of longer 
durations have been observed from several sources (e.g., 
Kaptein et al. (2000); in't Zand et al. (2002, 2005)). 
They have typical durations of « 30 minutes and ener- 
gies « 10 41 ergs. These bursts are known as long bursts. 
Even longer and more energetic thermonuclear bursts, viz. 
superbursts, have been observed from a few neutron star 
LMXBs (Cornelisse et al. (2000); Kuulkers (2004); in't 
Zand et al. (2004); Kuulkers (2005); Keck et al. (2008) 



may cause the superbursts (Cumming and Bildsten , 2001; 
Strohmayer and Brown , 2002). But this explanation has a 
problem for H/He accretors, because only a small amount 
of carbon is expected to remain after the burning of H and 
He via the rp-process. Although several solutions of this 
problem have been proposed by Cumming and Bildsten 

(2001) ; Cooper et al. (2006); Cooper et al. (2009), more 
detailed theoretical studies are required to understand the 
origin of superbursts. 

The study of thermonuclear bursts provides a unique 
opportunity to understand some aspects of extreme physics. 
In particular, these bursts can be very useful to constrain 
the neutron star parameters (see later sections) , and hence 
to understand the nature of the stellar core matter at 
supranuclear densities (see § 1). However, in order to use 
these bursts as a reliable tool, they have to be understood 
well (see, for example, Woosley et al. (2004); Weinberg 
et al. (2006a,b); Lin et al. (2006); Heger et al. (2007a); 
Cooper and Narayan (2007b); Fisker et al. (2008) for 
recent burst model calculations). One of the poorly ex- 
plored aspects of these bursts is the thermonuclear flame 
spreading on the neutron star surfaces. It is expected that 
the accreted matter quickly spreads all over the neutron 
star surface (Inogamov and Sunyaev , 1999). Therefore a 
burst may ignite at a certain location on the stellar surface 
(Cooper and Narayan , 2007a; Maurer and Watts , 2008), 
and then the burning region may expand to cover the en- 
tire surface (Fryxell and Woosley , 1982; Spitkovsky et al. 
, 2002). The study of this flame spreading is not only es- 
sential to understand the bursts, but also important to use 
these bursts as a tool to probe the neutron star parameters 
(§ 3.1.6). However, the theoretical modeling of thermonu- 
clear flame spreading including all the main physical effects 
is extremely difficult, and only recently Spitkovsky et al. 

(2002) have provided some insights in this field. Although 
these authors have ignored several physical effects (mag- 
netic field, strong gravity, etc.) in their model, they have 
considered the effects of the Coriolis force, which is im- 
portant as the bursting neutron stars are rapidly spinning 
(typical spin frequency <~ 300—600 Hz). The observational 
study of flame spreading is also no less difficult than the 
theoretical study. This is because, this spreading mostly 
happens within ss 1 sec during the burst rise, and in order 
to study it observationally, one needs to detect and mea- 
sure the evolution of the burst spectral and timing proper- 
ties during this short time. Therefore, although the obser- 
vational indications of flame spreading were found about a 
decade ago (e.g., Strohmayer et al. (1997b)), only recently 
such spreading could be confirmed and measured for a few 
bursts by analyzing the RXTE PCA data (Bhattacharyya 
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and Strohmayer , 2005, 2006a, b,c, 2007a, c). These findings 
provide the motivation for rigorous theoretical studies in- 
cluding all the main physical effects, as well as for detailed 
analysis of the burst rise data. 

3.1.2. Why are Thermonuclear X-ray Bursts Useful? 
Thermonuclear X-ray Bursts can be very useful to con- 
strain the neutron star parameters due to the following 
reasons. 

(1) They originate from the neutron star surfaces. There- 
fore, the stellar parameters (e.g., mass, radius, spin fre- 
quency) influence the burst spectral and timing proper- 
ties. Hence the modeling of these properties can be used 
to measure several neutron star parameters simultaneously 
(see later). 

(2) Various burst phenomena, such as photosphcric ra- 
dius expansion, thermonuclear flame spreading, etc., can 
be useful to constrain the neutron star parameters (see 
later). 

(3) Relatively low magnetic fields (<~ 10 7 — 10 9 G) of the 
bursting neutron stars simplify the modeling of the surface 
emission. 

(4) The surface emission during a burst is typically about 
10 times brighter than the rest of the X-ray emission. 
Hence the surface emission can be isolated from the to- 
tal emission relatively easily, and can be used to constrain 
the neutron star parameters with confidence. Note that 
this is not true for the non-burst emission, for which it is 
difficult to distinguish between the surface and the non- 
surface emissions with sufficient accuracy. 

(5) The ~ 10 times higher intensity of bursts, relative to 
the non-burst intensity, gives a significantly higher signal- 
to-noise ratio. 

(6) Many bursts can be observed from the same neutron 
star. The values of neutron star mass, radius and spin fre- 
quency, which affect the burst properties, do not change 
from one burst to another for a given source. Therefore, 
the joint analysis of many bursts, keeping the model values 
of these stellar parameters same for all the bursts, can be 
very useful to constrain these parameters. This technique 
is demonstrated in Bhattacharyya et al. (2005) (see also 
§ 3.1.8). 

(7) Thousands of thermonuclear X-ray Bursts have been 
observed from about 80 sources (Liu et al. (2007); see 
also Galloway et al. (2008a) for a catalog of RXTE PCA 
bursts). Therefore, a large amount of burst data is already 
available. 

3.1.3. Continuum Spectrum Method 

The continuum spectra of thermonuclear X-ray bursts 
are typically quite well described with a blackbody model 
(Swank et al. (1977); Hoffman et al. (1977b); Galloway 
et al. (2008a); see also Figure 3). The temperature of 
this blackbody increases rapidly as the accumulated mat- 
ter burns during the burst intensity rise, and then de- 
creases relatively slowly as the neutron star surface cools 
down during the burst decay. If the entire neutron star 



surface emits like a blackbody of the same temperature at 
a given time, then the inferred radius at infinity (Roo) of 
the neutron star can be obtained from the relation: 

Roo = (F^/aT^d. (2) 

Here is the observed bolometric flux, is the fit- 
ted blackbody temperature, d is the source distance and 
a is the Stcfan-Boltzmann constant. This provides one 
of the very few ways to measure the neutron star radius. 
Note that the radius measurement is typically more dif- 
ficult than the measurements of the neutron star mass 
and spin frequency. This is because, whereas the gravi- 
tational effect of mass influences the observed motion of 
the companion star in a binary system, and the observed 
flux may be modulated by the neutron star spin, the ra- 
dius exhibits no such effects. As a result, although the 
mass and the spin frequency of several neutron stars have 
been measured (Thorsett and Chakrabarty , 1999; Lamb 
and Boutloukos , 2008), accurate measurements of the ra- 
dius are challenging. 

The burst continuum spectrum method was used to 
constrain the neutron star radius soon after the discov- 
ery of these bursts (e.g., van Paradijs (1978); Goldman 
(1979); van Paradijs (1979, 1982); Marshall (1982); van 
Paradijs and Lewin (1986); Fujimoto and Taam (1986); 
Sztajno et al. (1987); Gottwald et al. (1987)). However, 
these efforts were not too successful because of the follow- 
ing systematic uncertainties. 

(1) Equation 2 assumes that the entire stellar surface emits 
during the burst. If that is not true, this equation will un- 
derestimate the value of the neutron star radius. 

(2) The lack of knowledge of the source distance introduces 
an uncertainty in the neutron star radius value inferred 
from this method. 

(3) The surface gravitational redshift (1 + z) makes the ob- 
served radius and the actual neutron star radius Rbb 
different by an unknown factor. 

(4) Although a blackbody model fits the observed burst 
spectrum well, this spectrum cannot be the blackbody 
originated from the burning layer. This is because the 
scattering of photons by the electrons and the frequency 
dependence of the opacity in the neutron star atmosphere 
hardens the energy spectrum (London et al. , 1984, 1986; 
Syunyaev and Titarchuk , 1986; Ebisuzaki and Nakamura , 
1988; Titarchuk , 1994; Madej et al. , 2004; Majczynaet al. 
, 2005). Such a hardening shifts the burning layer black- 
body spectrum towards the higher energies, and changes 
its shape slightly. This modified spectrum is called a "di- 
luted blackbody", and its temperature (T co i, i.e., the color 
temperature) is related to the actual blackbody tempera- 
ture of the burning layer (Tbb) by an unknown color factor 
/. Consequently, the observed radius is lower than the ac- 
tual radius by the square of the color factor. As a result, 
the actual temperature (Tbb) and radius (-Rbb) are related 
to the observed values by the following relations (includ- 
ing the effects of the gravitational redshift; Sztajno et al. 
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Figure 3: This figure shows the fitting of a thermonuclear burst spec- 
trum with an absorbed blackbody model. The spectrum is for a 1 sec 
segment of a burst observed from a neutron star LMXB 4U 1636-536 
with RXTE PCA. Since PCA cannot reliably measure the absorp- 
tion due to the Galactic neutral hydrogen because such absorption 
mostly happens outside the PCA band, we fixed the neutral hydro- 
gen column density njj to the estimated value (3.58 X 10 21 cm -2 ) 
in the source direction for the fitting. The upper panel shows the 
data points with error bars, and the best-fit model, while the lower 
panel shows the (data— model). This figure demonstrates that the 
continuum spectra of thermonuclear X-ray bursts are typically well 
described with a blackbody model (§ 3.1.3). 
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where the surface gravitational redshift 1+z = (1— 267 M/Rc 2 ' 
for a non-spinning neutron star. Here, R/M is the neutron 
star radius-to-mass ratio, G is the gravitational constant 
and c is the speed of light in vacuum. 

As mentioned earlier, thermonuclear bursts provide one 
of the very few promising methods to measure the neutron 
star radius. Therefore, although this method has a few sys- 
tematic uncertainties, it is worthwhile to try to improve it 
by removing or at least reducing these uncertainties. Here 
we briefly mention a few plausible ways to reduce the un- 
certainties. 

(1) Even when a part of the stellar surface burns, R 2 ^ 
should be proportional to, and hence is a measure of the 
burning region area. Hence, the uncertainty #1 should be 
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Figure 4: This figure shows the M — R space of neutron stars with the 
curves corresponding to a few representative EoS models (same as 
Figure 1). The green patch shows the allowed M — R space assuming 
Rbb = R and 1 + z = (1 - 2GM/Rc 2 ) -1 / 2 in equation 4, and for 
Roo-f 2 in the range 15 km — 18 km (see § 3.1.3 for details). This 
range should be reasonable, because Roc is often found in the 10 
km — 20 km range, and / can be within 1.2 and 1.8 according to 
Majczyna et al. (2005). This example figure demonstrates how the 
thermonuclear burst continuum spectra can be used to constrain the 
neutron star mass and radius, and hence the EoS models. 



reduced if the burning region area can be estimated inde- 
pendently; for example, by fitting the phase-folded burst 
oscillation light curves with a physical model (see § 3.1.6). 

(2) In order to reduce the uncertainty #2, one needs to 
measure the source distance by an independent method. 
This is possible if the source is in a globular cluster (e.g., 
Ozel et al. (2009)); because the distances of these clusters 
can typically be measured with sufficient accuracy. The 
source distance can also roughly be measured from the 
flux of a photospheric radius expansion burst (see § 3.1.5). 

(3) The surface gravitational redshift 1 + z depends on 
the stellar parameters. As mentioned earlier, for a non- 
spinning neutron star, 1 + z = (1 — 2GM/Rc 2 )~ 1 / 2 . There- 
fore, the uncertainty #3 can be reduced if the neutron 

_s i taj radius-to-mass ratio is constrained. Several obser- 
vational features of neutron star LMXBs provide inde- 
pendent methods to measure this ratio (for example, see 
§ 3.1.4; § 3.1.6; § 3.4.2). However, if the uncertainties 
#1, #2 and #4 can be sufficiently reduced, it is pos- 
sible to constrain the neutron star M — R space effec- 
tively, even without measuring the 1 + z. This is shown 
in Figure 4, and can be understood from equation 4 using 
l + z = (l-2G¥/iJc 2 )- 1 / 2 . 

(4) The color factor / depends on the neutron star surface 
gravity g, the actual blackbody temperature Tbb and the 
chemical composition in the stellar atmosphere (Madej et 
al. , 2004; Majczyna et al. , 2005). Therefore, in order 
to measure / and to reduce the uncertainty #4, one needs 
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to estimate these three parameter values, and model the 
atmosphere for these values. While the gravity g depends 
on the neutron star mass and radius (see equation 7; see 
§ 3.1.7 and Figure 10 for a plausible method to measure 
g), the latter two parameter values can be estimated from 
a realistic modeling of the thermonuclear bursts and the 
chemical composition of the burning layer. The determi- 
nation of the composition may be attempted using the 
following ways, (a) The observed duration of the burst, as 
well as the detailed modeling of its light curve and other 
properties may be useful to determine the chemical com- 
position (see § 3.1.1; also see Woosley et al. (2004); Heger 
et al. (2007a). (b) The observed nature of the companion 
star can be used to constrain the chemical composition of 
the accreted matter. Note that both the preburst burning 
layer and the preburst atmosphere consist of this matter, 
although the former may additionally contain ashes from 
the previous bursts. In case the companion is not directly 
observable, the measured binary orbital period may be use- 
ful to probe its nature. For example, if the system is an 
UCXB (orbital period < 80 min; in't Zand et al. (2008)), 
the companion and hence the accreted matter should be 
hydrogen deficient and heavier element rich (see § 3.1.1). 
(c) The observed atomic spectral lines from the non-burst 
emission can also be useful to estimate the composition 
of the accreted matter. Such lines are normally observed 
from high- inclination neutron star LMXBs (e.g., Jimenez- 
Garateetal. (2003); Diaz Trigo et al. (2006)). Therefore, 
the detailed observational and theoretical studies of burst 
and non-burst emissions can, in principle, constrain the 
burst color factor /. However, note that / is expected 
to change, as Tbb and plausibly the atmospheric chemical 
composition evolve during a burst, and hence one needs to 
estimate / vaules throughout the burst. 

The burning region is normally expected to cover the 
entire stellar surface during the burst decay (see § 3.1.1). 
If this is true, the uncertainty #1 does not appear, and 
the continuum spectrum of burst decay (rather than burst 
rise) should be more useful to measure the neutron star ra- 
dius. But unfortunately, the inferred burning area (oc i?^J 
can increase, decrease, or go up and down erratically in a 
given burst (see Figure 9 of Galloway et al. (2008a)). 
This may be because of the actual change of i? BB and/or 
the evolution of / (see equation 4) and/or any other rea- 
sons. Therefore, it is essential to understand the appar- 
ently erratic evolution of R 2 ^ for a reliable measurement 
of neutron star radii using burst decay continuum spec- 
tra. Very recently, Bhattacharyya et al. (2009) have, 
for the first time, reported a pattern in the apparently er- 
ratic behaviour of the inferred burning area. They found 
that the rate of change of the inferred area during decay 
is anticorrelated with the burst decay duration, and this 
anticorrelation is significant at the < 10~ 45 level (see Fig- 
ure 5). These authors have suggested that the variations 
in /, and more specifically the variations in the composi- 
tion of the atmosphere between bursts with long and short 
durations, may explain the anticorrelation. Such a varia- 
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Figure 5: Slope of burst blackbody radius-squared (-R^ oc inferred 
burning area) during burst decay vs. burst decay duration. Each 
point corresponds to one burst, and 877 bursts from 43 sources have 
been shown. This figure indicates a strong anticorrelation (significant 
at the < 10 — 45 level) between the rate of change of the inferred 
burning area and the decay duration, which can be useful to make 
the continuum burst spectrum method reliable (see § 3.1.3). For the 
exact definition of the radius-squared slope, see Bhattacharyya et al. 
(2009). 

tion can be reasonable, because the duration of a burst 
depends on the chemical composition of the burning ma- 
terial (see § 3.1.1). However, in order to verify this model, 
extensive theoretical study is required. But, whatever is 
the correct model, the newly discovered pattern in the R 2 ^ 
evolution will play an extremely important role to make 
the continuum spectrum method reliable. 

3.1.4- Spectral Line Method 

Observation and identification of narrow atomic spec- 
tral lines from the surface of a neutron star (for exam- 
ple, during a thermonuclear burst) provide the cleanest 
way to measure the stellar radius-to-mass ratio using the 
following formula (for the Schwartzschild spacetime ap- 
propriate for a non-spinning neutron star): E cm /E b s = 
1 + z = [1 - {2GM / Re 2 )]- 1 / 2 (Ozel and Psaltis (2003); 
Bhattacharyya et al. (2006a); see Figure 6). Here, E cm 
is the emitted energy of the line photons, and E b s is 
the observed energy of these photons. However, if the 
surface lines are broad and asymmetric, E b s cannot be 
determined uniquely, and hence it is difficult to measure 
Rc 2 /GM accurately using the above formula. Unfortu- 
nately, neutron stars in LMXB systems typically spin with 
high frequencies (~ 300 — 600 Hz), and hence the Doppler 
effect and the special relativistic beaming are expected to 
make the surface lines broad and asymmetric. Recently, 
Bhattacharyya et al. (2006a) have suggested a way to 
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measure the stellar Rc 2 /GM using the above formula for 
such broad lines. They have shown that using y/E\Eh as 
E'obs, one can measure Rc 2 /GM with less than 2% error. 
Here, E\ is the energy of the low energy end of the broad 
line, and E^ is that of the high energy end. This shows 
that surface lines can be used to accurately measure the 
stellar Rc 2 /GM even for LMXBs having rapidly spinning 
neutron stars. 



easier by keeping the magnetic splitting negligible; and (4) 
the luminous bursts give good signal-to-noise ratio. More- 
over, the modeling of the surface lines detected from ther- 
monuclear bursts may be easier, because (1) magnetic field 
is relatively low; and (2) during the bursts, typically 90% 
of the total emission originates from the neutron star sur- 
face, and hence the surface spectral line should be largely 
free from the uncertainty due to the other X-ray emission 
components, such as the accretion disk. 




Radius (km) 



Figure 6: This figure shows the M — R space of neutron stars with 
the curves corresponding to a few representative EoS models (same 
as Figure 1). The green patch shows the allowed M — R space for a 
reasonable range of the radius-to-mass ratio {Rc 2 /GM = 4.0 - 4.2). 
This example figure demonstrates how the measured radius-to-mass 
ratio can be used to constrain the neutron star M — R space (see 
§ 3.1.4; § 3.1.6; and § 3.4.2 for details). 

A surface absorption line originated from the entire 
stellar surface can have only one dip (or peak for an emis- 
sion line; see Figure 1 of Ozel and Psaltis (2003)). On the 
other hand, if it originates from a portion of the surface, 
especially near the equator, it can have two dips (or two 
peaks for an emission line; Ozel and Psaltis (2003); Bhat- 
tacharyya et al. (2006a); see also Figure 7). In either 
case, its width and shape depends on the neutron star pa- 
rameter values (Figure 7; see also Ozel and Psaltis (2003); 
Bhattacharyya et al. (2006a)). Therefore, fitting the sur- 
face lines with physical models can be useful to constrain 
additional (apart from Rc 2 /GM) neutron star parameters. 

Surface lines are likely to be detected and identified 
from thermonuclear X-ray bursts and hence from neutron 
star LMXB systems. This is because: (1) the continuous 
accretion and the radiative pressure may keep the line- 
forming heavy elements in the stellar atmosphere for the 
time-duration required for the line detection; (2) nuclear 
burning ashes may be transported to the atmosphere by 
strong convection (Weinberg et al. , 2006a); (3) the rela- 
tively low magnetic fields (~ 10 7 — 10 9 Gauss) of the neu- 
tron stars in LMXBs should make the line identification 



l.oo r 

g 0.98 - 




0.60 0.65 0.70 0.75 0.00 0.85 0.90 0.95 

E ota /E em 



Figure 7: Profiles of narrow absorption lines formed at the equa- 
torial portion of the neutron star surface as they would appear to 
an observer at i sp i n = 75°, where i ap i n is the observer's inclination 
angle measured from the stellar spin axis (see Bhattacharyya et al. 
(2006a)). Here, E e m is the emitted energy, -E bs ls the observed en- 
ergy, and the continuum is assumed to be at a constant photon flux 
of 1 (in arbitrary unit). The solid profile is for the neutron star pa- 
rameter values M = IAMq and R = 9 km; and the dotted profile is 
for M = 1.5Mq and R = 11 km. All other source parameter values 
are same, including the stellar spin frequency (= 400 Hz). This fig- 
ure shows that the shape and the width of the surface spectral lines 
depend on the neutron star parameters (§ 3.1.4). 

In fact, the only plausible observation of the neutron 
star surface atomic spectral lines was from an LMXB EXO 
0748-676 (Cottam et al. , 2002). These authors detected 
two significant absorption features in the XMM-Newton 
Reflection Grating Spectrometers (RGS) energy spectra of 
the thermonuclear X-ray bursts, which they identified as 
Fe XXVI and Fe XXV n = 2 — 3 lines with a gravitational 
redshift of 1+z = 1.35. In order to detect these lines, they 
had to combine the spectra of 28 bursts. The absorption 
of some of the burst photons (originated in deep burning 
layers) by the iron ions in the upper atmosphere might 
give rise to these redshifted lines. Chang et al. (2005) 
showed that the observed strength of the iron lines could 
be produced by a neutron star photospheric metallicity, 
which was 2 — 3 times larger than the solar metallicity (see 
also Chang et al. (2006)). The identification of these 
two lines as the surface atomic spectral lines is reason- 
able for the following reasons: (1) the early burst phases 
(when the temperature was higher) showed the Fe XXVI 
line, while the late burst phases (when the temperature 
was lower) exhibited the Fe XXV line, which is consistent 
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with the qualitative temperature dependence of the iron 
ionization states; (2) the narrow lines were consistent with 
the relatively slow stellar spin rate (45 Hz), reported by 
Villarreal and Strohmayer (2004) (but see later); and (3) 
both the lines showed the same surface gravitational red- 
shift, and this redshift value is reasonable for the surface of 
a typical neutron star. However, recent observations have 
argued against the neutron star surface origin of these two 
lines. (1) these lines were not significantly detected from 
the later observations of EXO 0748-676 (Cottam et al. , 
2008). However, such disappearence might be caused by 
the change in the photospheric conditions, which might 
weaken the lines. (2) Galloway et al. (2009) have discov- 
ered burst oscillations at 552 Hz from EXO 0748-676 with 
a significance of 6.2 a. This argues that the neutron star 
spin frequency is 552 Hz (see § 3.1.6). If this is true, the 
observed narrow spectral lines from EXO 0748-676 could 
originate from the stellar surface only if the entire line 
forming region were very close to the neutron star spin 
axis (Bhattacharyya et al. , 2006a). Moreover, recently 
Balman (2009) has suggested that the 45 Hz frequency re- 
ported by Villarreal and Strohmayer (2004) may be a vari- 
ability within the boundary layer. Therefore, the surface 
origin of the observed absorption features from EXO 0748- 
676 is not established yet. Besides, Kong et al. (2007) did 
not find surface spectral lines from the combined XMM- 
Newton RGS spectra of 16 thermonuclear bursts from the 
neutron star LMXB GS 1826-24. 

3.1.5. Photospheric Radius Expansion Burst Method 

As mentioned in § 3.1.1, for some of the strong bursts, 
viz. photospheric radius expansion (PRE) bursts, the pho- 
tospheric radiative pressure force exceeds the gravitational 
force, and the photosphere is temporarily lifted off the 
neutron star surface. During the expansion of the photo- 
sphere, the temperature decreases and the burning area 
inferred from the continuum spectroscopy increases. Sub- 
sequently, as the photosphere comes down, the temper- 
ature increases and the inferred burning area decreases. 
After the photosphere settles on the stellar surface at the 
touchdown time, it starts cooling down, and the temper- 
ature decreases again. As a result, PRE bursts can be 
identified by two maxima in the temperature profile and 
a burning area maximum corresponding to the tempera- 
ture minimum (see Figure 9 of Galloway et al. (2008a)). 
The amount of photospheric expansion can vary by a large 
extent from burst to burst. Normally during the expan- 
sion, the X-ray flux does not change much. But for cer- 
tain very powerful bursts the expansion can be so large 
that the photospheric temperature may be shifted below 
the X-ray band, and the burst may be divided into two 
parts: a precursor and a main part (Lewin et al. , 1984; 
Strohmayer and Brown , 2002). For other somewhat less 
powerful PRE bursts, a part of the flux may shift out of 
the X-ray band making the burst intensity profile double- 
peaked (Strohmayer and Bildsten , 2006). 

The Eddington limit (see § 3.1.1) on the neutron star 
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Figure 8: This figure shows the M — R space of neutron stars with 
the curves corresponding to a few representative EoS models (same 
as Figure 1). The green patch shows the allowed M — R space using 
equation 6 (see § 3.1.5 for details). In this example figure we have en- 
sured that the green patch passes through the the reasonable ranges 
of M and R, and this figure demonstrates how the photospheric ra- 
dius expansion bursts can be used to constrain the neutron star mass 
and radius, and hence the EoS models. 

surface is given by 

iEdd = (4ttGMc/k)(1 - 2GM/Rc 2 y 1/2 = Ai:R 2 aT^ B , (5) 

where, k is the atmospheric opacity (e.g., Strohmayer and 
Bildsten (2006)). Here we have assumed that the entire 
stellar surface emits, and the surface gravitational redshift 
1 + z = [1 — (2GM/Rc 2 )}^ 1 ^ 2 , which is appropriate for a 
non-spinning neutron star. It may be reasonable to iden- 
tify the observed touchdown flux as the Eddington-limited 
flux (Ozel , 2006). Therefore, equations 5 and 3 for the 
touchdown time can be used to constrain the neutron star 
M — R space, if the entire stellar surface emits (see Gal- 
loway et al. (2003)) and the parameters k and color factor 
/ are known. The flux observed at infinity corresponding 
to equation 5 is (e.g., Ozel (2006)) 

^Edd = (l/4vrd 2 )(47rGMc/K)(l - 2GM/Rc 2 ) 1/2 , (6) 

where d is the source distance. Considering the observed 
touchdown flux as i<Edd (Ozel , 2006), equation 6 can be 
used to constrain the neutron star M — R space (see Fig- 
ure 8), if the entire stellar surface emits (Galloway et al. 
, 2003) and the parameters d and k are known. In order 
to estimate k and/or /, one needs to model the stellar at- 
mosphere for a known chemical composition. The possible 
ways to determine this composition is discussed in § 3.1.3. 
Apart from the above methods, the observed variation of 
the Eddington luminosity and the corresponding change in 
the gravitational redshift caused by the photospheric ex- 
pansion may, in principle, be useful to determine the neu- 
tron star surface gravitational redshift (sec for example, 
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Damen et al. (1990); van Paradijs et al. (1990)). These 
and related methods have been used by many scientists in 
order to constrain the neutron star mass and radius (e.g., 
Fujimoto and Taam (1986); Ebisuzaki (1987); Sztajno et 
al. (1987); van Paradijs and Lewin (1987); Chevalier and 
Ilovaisky (1990); Kaminkcr et al. (1990); Damen et al. 
(1990); van Paradijs et al. (1990); Haberl and Titarchuk 
(1995); Smale (2001); Strohmayer et al. (1998); Titarchuk 
and Shaposhnikov (2002); Shaposhnikov et al. (2003); 
Shaposhnikov and Titarchuk (2004); Ozcl (2006)). How- 
ever, so far such measurements heve not been generally 
precise enough, because of the systematic uncertainties 
primarily due to the unknown parameter values mentioned 
above (e.g., Strohmayer and Bildstcn (2006)). 

This and other sections (e.g., § 3.1.3; § 3.5.2) point out 
that the knowledge of the source distance d can be very 
useful to measure the neutron star parameters. In 1978, 
van Paradijs (1978) suggested that the photospheric ex- 
pansion flux, which corresponds to the Eddington flux, can 
be a "standard candle", and hence a distance indicator 
(e.g., using equation 6). In the next 30 years, scientists 
have tested this idea by examining the variation of the 
photospheric expansion flux from burst to burst, compar- 
ing the inferred distance to the known distance (whenever 
available; e.g., for globular cluster sources), and in other 
ways (e.g., van Paradijs (1981); Basinska et al. (1984); 
Smale (1998); Kuulkers and van der Klis (2000); Kuulkers 
etal. (2003); Galloway et al. (2003, 2006); Galloway and 
Cumming (2006); Ozcl (2006); Galloway et al. (2008b); 
Table 9 of Galloway et al. (2008a)). They have found 
that this method is affected by the uncertainties due to 
unknown neutron star mass, surface gravitational redshift, 
photospheric composition, etc., and hence the distance can 
typically be constrained within a somewhat large range. 

3.1.6. Burst Oscillation Method 

High frequency narrow timing features, viz. burst oscil- 
lations, have been detected from many thermonuclear X- 
ray bursts with the PCA instrument of RXTE (Strohmayer 
and Bildsten , 2006). They can be observed either in the 
burst rising part or in the decay part, or in both. How- 
ever, neither these oscillations appear in every burst from 
a given source, nor they have been observed from every 
neutron star LMXB (Galloway et al. , 2008a). Burst os- 
cillations were first discovered from the LMXB 4U 1728-34 
by Strohmayer et al. (1996) in 1996. After that, this fea- 
ture has been detected from 20 other neutron star LMXBs 
(oscillations from three of them are not yet confirmed; see 
Bhattacharyya (2007); Lamb and Boutloukos (2008); Gal- 
loway et al. (2009); Watts et al. (2009a); Markwardt et 
al. (2007); Watts et al. (2009b) and references therein). 

Burst oscillations can be seen as a peak in the Fourier 
power spectrum computed from a high time resolution 
light curve of either the entire burst or a part of it (see 
Figure 9). The frequency of this feature can evolve during 
a burst (Strohmayer and Bildstcn (2006) and references 
therein). Normally frequency shifts by less than 1% of the 
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Figure 9: Upper panel: X-ray intensity profile of a thermonuclear X- 
ray burst from the neutron star LMXB 1A 1744-361 (Bhattacharyya 
et al. , 2006b). Power contours around 530 Hz implies that there are 
significant intensity variations (burst oscillations; § 3.1.6) with as 530 
Hz frequency. Note that the shape of the power contours indicates 
a slight change in frequency with time, which is common for burst 
oscillations. Lower panel: The corresponding power spectrum. The 
peak at fs 530 Hz implies the burst oscillations with that frequency. 
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mean frequency, but it can be somewhat greater for some 
bursts. The frequency normally increases from a lower 
initial value to an asymptotic value, but it has also been 
seen to decrease for some bursts (see for example, Muno et 
al. (2002a)). It is believed that burst oscillations originate 
from an azimuthally asymmetric brightness pattern on the 
spinning neutron star surface. Since in such a case the X- 
ray flux measured by the observer should vary with the 
stellar spin frequency, the burst oscillation (asymptotic) 
frequency is believed to be the stellar spin frequency. This 
model is established because of the following reasons. 

(1) Burst oscillations originate from the thermonuclear 
bursts which are neutron star surface phenomena. 

(2) The oscillation (asymptotic) frequency is very stable 
on timescales of years for a given source (Muno ct al. , 
2002a). This strongly suggests that burst oscillations are 
connected to the stellar spin, which is also expected to be 
a stable property. 

(3) The observed burst oscillation frequencies are mostly 
in the range ~ 300 — 600 Hz. Such rapid spin rates are ex- 
pected for neutron stars in LMXBs, because of the accretion- 
induced angular momentum transfer (§ 1). 

(4) Spin frequencies of five burst oscillation sources are 
independently known from accretion-powered pulsations 
(see § 3.2). For each of them the burst oscillation fre- 
quency matches the known frequency well (Chakrabarty 
et al. , 2003; Strohmayer et al. , 2003; Casella et al. , 
2008; Watts et al. , 2009a,b). 

The evolution of burst oscillation frequency may be caused 
by the slow motion of the photosphere on the stellar sur- 
face. One reason for such a slow motion may be the angu- 
lar momentum conservation of contracting thermonuclear 
shells (Strohmayer et al. , 1997a). In this scenario, a ther- 
mally expanded shell has a spin frequency lower than the 
stellar spin frequency due to angular momentum conser- 
vation. As the shell contracts, its spin frequency increases 
and tends to the stellar spin frequency. However, this 
model cannot explain a frequency shift roughly larger than 
2 Hz (Strohmayer and Bildsten , 2006) observed from sev- 
eral bursts (e.g., Galloway et al. (2001); Wijnands et al. 
(2001); Chakrabarty et al. (2003)). Therefore, other phys- 
ical reasons (e.g., involving magnetic field; Bhattacharyya 
and Strohmayer (2006c); Heng and Spitkovsky (2009)) 
may contribute to the frequency evolution process. 

Burst oscillations provide one of the two ways to mea- 
sure the spin rates of neutron stars in LMXB systems (see 
§ 3.2 for the other method). In fact, this feature provided 
most of the known frequencies of such neutron stars. In 
addition, the amplitude and the shape of a phase-folded 
burst oscillation light curve can be very useful to constrain 
other stellar parameters. This is because the compactness 
(GM/Rc 2 ) and the spin- related surface speed (i?) affect 
these light curves through the Doppler effect, and special 
and general relativistic effects (including light bending). 
The amplitude and asymmetry of phase- folded light curves 
are greater for lesser GM/Rc 2 values (e.g., Pechenick et 
al. (1983); Strohmayer et al. (1998); Miller and Lamb 



(1998); Brajeet al. (2000); Weinberg et al. (2001)), and 
hence can be used to constrain the neutron star radius- 
to-mass ratio. Besides, i9 introduces a pulse phase depen- 
dent Doppler shift in the X-ray spectrum (Strohmayer and 
Bildsten , 2006), which can be used to estimate d. Note 
that $ can be useful to constrain the neutron star radius, 
for the known stellar spin frequency (measured from burst 
oscillations). 

Several studies have been done to measure the neutron 
star properties based on the modeling of burst oscillations. 
These include an investigation of the amplitude and har- 
monic content of the phase-folded light curves due to a 
point-like hot spot (Miller and Lamb , 1998); and the ef- 
fects of expanding single hot spot or two hot spots (Nath 
et al. , 2002), or single and antipodal hot spots of varying 
size (Weinberg et al. , 2001). Moreover, Cadeau et al. 
(2007) and Morsink ct al. (2007) studied the effects of 
the oblate shape of a rapidly spinning neutron star on the 
phase-folded light curves. 

Apart from XTE J1814-338, no other neutron star LMXB 
has exhibited significantly asymmetric (nonsinusoidal) phase- 
folded burst oscillation light curves (but see Bhattacharyya 
and Strohmayer (2005)). The symmetric (sinusoidal) na- 
ture of the phase-folded light curves can be attributed to 
the limitation of the detector capability and/or to one of 
the following reasons (Muno et al. , 2002b): (1) if there is 
one hot spot on the neutron star, it must either be close 
to the spinning pole or cover nearly half the stellar sur- 
face; or (2) if there are two hot spots, they must form 
very near the spinning equator. Information about the 
neutron star parameters can be extracted only from the 
amplitude of such a symmetric light curve, and not from 
its shape. Therefore, Bhattacharyya et al. (2005) fitted 
the nonsinusoidal burst decay oscillation light curves of 
XTE J1814-338 with a relativistic model appropriate for a 
single hot spot. This model had five input parameters for 
a chosen EoS model and a known neutron star spin fre- 
quency. These parameters were, (1) neutron star radius- 
to-mass ratio (Rc 2 /GM), (2) polar angle of the center of 
the spot, (3) angular radius of the spot, (4) a parame- 
ter which gave a measure of the beaming in the emitter's 
frame, and (5) the observer's inclination angle measured 
from the stellar spin axis. Bhattacharyya et al. (2005) 
could determine a lower limit of Rc 2 /GM at the 90% confi- 
dence level, in spite of the systematics due to the unknown 
values of the other parameters. Watts et al. (2005) and 
Watts and Strohmayer (2006) studied the variability and 
the energy dependence of the burst oscillation properties 
of XTE J1814-338, which may be useful to reduce these 
systematic uncertainties. 

Since XTE J1814-338 is an accretion-powered millisec- 
ond pulsar (§ 3.2.1), and its modulations in the accretion 
and thermonuclear burst emission are phase-locked with 
each other (see Watts et al. (2008)), its burst decay os- 
cillation could very well originate from a hot spot. But 
it is difficult to understand what gives rise to a hot spot 
for a non-pulsar during burst decay. Therefore, several 
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groups have attempted to model burst decay oscillations of 
non-pulsars using an asymmetric brightness pattern more 
complex than a hot spot. Such a complex pattern may be 
caused by Rossby waves (or r-modes) in the surface layers 
(Lee and Strohmayer , 2005; Heyl , 2005; Piro and Bildsten 
, 2006), shear instabilities (Cumming , 2005), etc. Neutron 
star parameters (e.g., Rc 2 /GM) can, in principle, be con- 
strained, even if burst oscillations originate from surface 
modes (see Lee and Strohmayer (2005); Heyl (2005)). 
However, in practice, it is very difficult, because one needs 
to identify the burst oscillations with the correct surface 
mode number. 

Unlike burst decay oscillations of non-pulsars, the os- 
cillations during burst rise are likely to be caused by a hot 
spot (e.g., Bhattacharyya and Strohmayer (2007c); Watts 
et al. (2009a)). This hot spot is an expanding burn- 
ing region, which can be modelled more reliably than the 
burst decay oscillations, in order to constrain the neutron 
star parameters. A better understanding of thermonuclear 
flame spreading (§ 3.1.1) will further increase the reliabil- 
ity of neutron star parameter measurement using burst 
rise oscillations. However note that, although an expand- 
ing hot spot is likely to cause the burst rise oscillations, 
this hypothesis is yet to be comprehensively tested against 
the entire data set. 

3.1.7. Millihertz Quasi-periodic Oscillation Method 

In 2001, Revnivtsev et al. (2001) discovered milli- 
Hertz (mHz) period quasi-periodic oscillations (QPOs) from 
three neutron star LMXBs: 4U 1608-52, 4U 1636-536 and 
Aql X-l. This new class of timing feature has the following 
properties: (1) the frequency varies between 0.007— 0.009 
Hz (period « 1.9 — 2.4 minutes); (2) the associated flux 
variations are at the few percent level; (3) the strength is 
more at lower photon energies (< 5 keV; this is in contrast 
to other types of QPOs); (4) it appears in a narrow lumi- 
nosity range around 10 37 ergs s _1 ; and (5) it disappears 
after thermonuclear X-ray bursts. Based on the theoretical 
research and observational findings of Paczynski (1983); 
Bildsten (1993); Revnivtsev et al. (2001); Yu and van der 
Klis (2002); Narayan and Heyl (2003), and others, Heger 
et al. (2007b) proposed that the mHz QPOs correspond to 
the marginal thermonuclear burning stability on the neu- 
tron star surfaces close to the bounday between regime 3 
and regime 4 mentioned in § 3.1.1. Near such a bounday, 
oscillations are expected because the eigenvalues of the 
system are complex (Paczynski , 1983; Narayan and Heyl 
, 2003). The narrow range of accretion rate corresponding 
to the bounday can explain why mHz QPOs appear in a 
narrow luminosity range. Heger et al. (2007b) found that 
the period of the oscillations due to the marginal burning 
stability should be « v^thcrmai^accrction (~ 2 min), where 
^thermal is the thermal timescale of the burning layer and 
^accretion is the time to accrete the fuel at the Eddington 
accretion rate. This theoretical period matches well with 
the observed period, which is the main strength of the 
Heger et al. (2007b) model. However, note that a few 



aspects of mHz QPOs are still not understood (Heger et 
al. , 2007b). 
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Figure 10: This figure shows the M — R space of neutron stars with 
the curves corresponding to a few representative EoS models (same 
as Figure 1). The green patch shows the allowed M — R space using 
equation 7 for a reasonable range of the surface gravity (g = [2.6 — 
3.1] X 10 14 cm s — 2 ). This example figure demonstrates how the 
measured surface gravity can be used to constrain the neutron star 
M — R space (see § 3.1.7). 

The theoretically proposed oscillation frequency is a 
strong function of the surface gravity and the chemical 
composition of the accreted matter (see Figure 9 of Heger 
et al. (2007b)). Therefore, if the Heger et al. (2007b) 
model is correct, and if the chemical composition can be 
determined (see § 3.1.3 for a discussion of plausible meth- 
ods), then mHz QPOs can be used to constrain the M — R 
space of neutron stars (see Figure 10). This is because the 
surface gravity g is a function of stellar mass and radius 

g ~ GM/R 2 [\ — 2GM/Rc 2 ]~ 1 / 2 . (7) 

Here we note that the marginal burning stability depends 
on the stellar spin rate (Keek et al. , 2009), and hence the 
neutron star spin can also affect the mHz QPOs. However, 
fortunately the stellar spin frequencies of all the currently 
known mHz QPO sources have meen measured (Lamb and 
Boutloukos , 2008). 

Yu and van der Klis (2002) found that kilohertz QPO 
(see § 3.3) frequency is anticorrelated with the luminosity 
when mHz QPOs are observed. This is opposite to the 
long-term trend. Therefore, mHz QPOs affect the kilo- 
hertz QPOs, and hence can be useful to constrain the mod- 
els of the latter. Such a constraint can be very helpful to 
measure the neutron star parameters (see § 3.3.2; § 3.4.2). 
Altamirano et al. (2008) found that mHz QPOs can also 
be used to predict the occurrence of thermonuclear bursts. 
Such a prediction may be helpful to observe more bursts, 
as the duty cycle of bursts is small. 
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3.1.8. Joint Application of the Methods based on Ther- 
monuclear X-ray Bursts 

Thermonuclear X-ray bursts provide various methods 
to measure the neutron star parameters (§ 3.1.3; § 3.1.4; 
§ 3.1.5; § 3.1.6; § 3.1.7). Figures 4, 6 and 8 show how 
some of these methods can be used to constrain the neu- 
tron star M — R space. Several of the burst properties, 
which provide these methods, can be observed from a 
given neutron star LMXB. For example, apart from the 
surface spectral lines, all these burst features are observed 
from 4U 1636-536. Similarly, all the burst features, ex- 
cept mHz QPOs, have possibly been observed from EXO 
0748-676. Therefore, joint application of the burst meth- 
ods can be very useful to constrain the parameters of the 
same neutron star. For example, the continuum spectrum 
method (§ 3.1.3), the spectral line method (§ 3.1.4) and 
the PRE burst method (§ 3.1.5) together may constrain 
the M — R space effectively (e.g., the intersection of the 
green patches of Figures 4, 6 and 8; see also Figure 1 
of Ozel (2006)). The first and the third of the above 
methods, and the known source distance, can also reduce 
the systematic uncertainties (e.g., Ozel et al. (2009)). 
Moreover, Bhattacharyya and Strohmayer (2007a) have 
demonstrated that the joint application of the continuum 
spectrum method and the burst oscillation method (§ 3.1.6) 
can constrain the neutron star parameters by eliminating 
the source distance parameter. 

A given neutron star LMXB can exhibit many bursts 
with various properties (e.g., oscillations, photospheric ra- 
dius expansion, etc.). These bursts may be observed with 
different satellites (e.g., RXTE: XMM-Newton) at differ- 
ent times. The burst properties depend on neutron star 
and other system parameters. Neutron star global param- 
eter values remain unaltered for all these observations, al- 
though some other source parameter values (e.g., the hot 
spot size for oscillations) are expected to change. Here 
we briefly describe how these observations can be com- 
bined within the framework of Bayesian statistics, in order 
to constrain the neutron star parameters effectively (see 
Bhattacharyya et al. (2005)). Let the posterior probabil- 
ity density (which we determine by fitting the data with 
a model) over the full set of model parameters x\...x n be 
p{x\...x n ). If we are interested in the confidence regions 
for a single parameter Xk, then we integrate this probabil- 
ity distribution over the "nuisance" parameters x\...Xk-i 
and Xk+i-.x n : 

q{x k ) = J dx 1 ...dxk-idxk+i...dx n p(x 1 ...x n ). (8) 

The maximum likelihood of as well as its confidence re- 
gions, can be obtained from the probability function q(xk). 
Using this method, one can obtain a tight joint constraint 
on a neutron star global parameter (e.g., mass, radius, 
etc.) which does not vary from one observation to an- 
other, by integrating over the nuisance parameters of all 
the observations. 



3.2. Accretion-powered Millisecond Period Pulsations 
3.2.1. What is an Accretion-powered Millisecond Pulsar? 

Several neutron star LMXBs exhibit coherent fast (pe- 
riod ~ 1 millisecond) brightness pulsations in their persis- 
tent X-ray emissions (e.g., non-burst, non-quiescent, etc.; 
§ 3.1.1; § 3.5.1; wijnands (2005); Poutanen (2008); Al- 
tamirano et al. (2009a) and references therein). These 
pulsations can be detected as a narrow peak in the power 
spectrum. In these sources, the accretion flow is channeled 
by the neutron star magnetic field to the magnetic poles. 
The X-ray hot spots on the neutron star created by this 
channeled flow rotate around the stellar spin axis (if the 
magnetic axis and the spin axis are misaligned) and give 
rise to the observed X-ray pulses. Note that one of the 
hot spots is usually obsecured by the accretion disk. Since 
the gravitational energy release from the accreted mat- 
ter powers the X-ray pulsations, these sources are called 
"accretion-powered millisecond pulsars (AMPs)". 

In 1998, the first AMP SAX J1808.4-3658 was discov- 
ered with the RXTE satellite (Wijnands and van der Klis 
, 1998a; Chakrabarty and Morgan , 1998). This AMP 
exhibits significant pulsations at w 401 Hz. The data re- 
vealed that (1) the low-energy pulses lag the high-energy 
pulses, and (2) the pulse profiles are significantly asymmet- 
ric (nonsinusoidal) allowing a detection of even the second 
harmonic (Cui et al. , 1998). Poutanen and Gierlihski 
(2003) suggested that both these properties could be ex- 
plained in terms of two spectral components: a 'pencil'-like 
beamed blackbody emission and a 'fan'-like Comptonized 
emission (from a radiative shock; but see also Hartman et 
al. (2009)). 

Four years after the discovery of the first AMP, the 
second AMP XTE J1751-305 was reported by Markwardt 
et al. (2002). Since then ten more AMPs have been dis- 
covered (wijnands (2005); Altamirano et al. (2009a, b); 
Markwardt et al (2009); and references therein). All these 
AMPs, including the first two, are transient sources (pul- 
sations are detected only during the outbursts; § 3.5.1; 
Lamb et al. (2008a)), and three of them show pulsations 
very occasionally (see Watts et al. (2009a) and references 
therein) . Here we do not describe the properties of these 
AMPs in detail, and suggest the following references for 
the interested readers: Gilfanov et al. (1998); Psaltis and 
Chakrabarty (1999); Bildsten and Chakrabarty (2001); 
Miller et al. (2003); Nelson and Rappaport (2003); wi- 
jnands (2005); van Straaten et al. (2005); Poutanen 
and Beloborodov (2006); van der Klis (2006); Watts et 
al. (2008); Poutanen (2008); Cackett et al. (2009a); 
etc. Note that, although many properties of AMPs are 
broadly understood, it is still a matter of debate why only 
a few neutron star LMXBs arc AMPs (e.g., Titarchuk et al. 
(2002); Lamb et al. (2008a); Ozel (2009)). Study of the 
intermittent pulsation sources, among other things, may 
be useful to solve this problem (e.g., Lamb et al. (2008b); 
Watts et al. (2009a)). 
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3.2.2. Accretion-powered Millisecond Pulsation Method 

Accretion-powered pulsation frequency is the neutron 
star spin frequency, because such pulsations are caused 
by the periodic motion of a hot spot due to the stellar 
spin (see 3.2.1). These pulsations are more coherent than 
burst oscillations (§ 3.1.6), and, unlike the burst oscillation 
frequency, the pulsation frequency does not evolve appre- 
ciably in a short time. Therefore, the pulsation method 
gives a more accurate measurement of the spin frequency. 

The phase- folded accretion-powered pulsation light curves 
are affected by the Doppler effect, and the special and gen- 
eral relativistic effects. Therefore, similar to the burst os- 
cillations, the amplitude and the shape of a pulsation light 
curve can also be very useful to measure the neutron star 
and other system papameters (see § 3.1.6 for the burst os- 
cillation method). Various groups attempted to constrain 
the stellar parameters by fitting the pulsation light curves 
with appropriate relativistic models (e.g., Poutanen and 
Gierlihski (2003, 2004); Leahy et al. (2008, 2009)). How- 
ever, similar to the burst oscillation method (see § 3.1.6), 
the systematic uncertainties due to the unknown system 
parameters (such as, observer's inclination angle, size and 
polar angle of the hot spot, etc.) hinder an accurate mea- 
surement. Besides, additional systematics due to the un- 
certainties in the spectral properties of various X-ray emit- 
ting components make the pulsation method less reliable 
(see the next paragraph). 

Both burst oscillations and accretion-powered pulsa- 
tions originate from azimuthally asymmetric brightness 
patterns on the spinning neutron star surfaces. There- 
fore, in principle, similar methods can be used to extract 
information from them. However, in practice, these two 
methods are somewhat different from each other because 
of the following reasons. 

(1) Phase- folded burst oscillation light curves appear sym- 
metric or sinusoidal with the current observational capa- 
bilities (see, Strohmayer et al. (2003) and Bhattacharyya 
and Strohmayer (2005) for exceptions). However, the 
accretion-powered pulsation light curves can be nonsinu- 
soidal, especially, but not exclusively, at higher energies 
(e.g., Cui et al. (1998); Poutanen (2008)). Therefore, 
in principle, more information can be extracted from the 
accretion-powered pulsation light curves (using both am- 
plitudes and shapes; see § 3.1.6). 

(2) During a thermonuclear burst, most of the emission 
originates from the neutron star surface, and the corre- 
sponding energy spectrum can be fitted well with a black- 
body model (sec § 3.1.3). Moreover, only burst emission 
contributes to the oscillations, and the gravitational en- 
ergy release due to the accretion can be separated from 
this burst emission (see Galloway et al. (2008a)). These 
make the modeling of phase-folded burst oscillation light 
curves simple and reliable. But during accretion-powered 
pulsations, gravitational potential energy powers both the 
disk and the surface emissions. Moreover, the surface emis- 
sion, which gives rise to the pulsations, is at least partially 



processed (Comptonized) in the bulk motion of matter 
(e.g., Mastichiadis and Kylafis (1992)) channelled by the 
stellar magnetic field, and/or in a radiative shock at the 
neutron star surface (Poutanen and Gierlihski , 2003). It 
is very difficult to observationally separate these various 
components of emission, as well as their contributions to 
the phase-folded pulsation light curves, from each other. 
This makes the understanding of pulsation light curves dif- 
ficult, and reduces the reliability of the pulsation method. 
Rigorous theoretical computation of both time-averaged 
and phase-resolved energy spectra is required to solve this 
problem. For example, among other things, one may need 
to perform the Monte Carlo simulation of the scattering of 
photons in the bulk electron flow near the magnetic poles, 
including the gravitational redshift and light bending ef- 
fects. 

(3) Unlike burst oscillations, accretion-powered pulsations 
can be useful to test the effects of accretion rate changes 
on the inner-edge radius of the accretion disk, the size of 
the hot spot, etc. (Ibragimov and Poutanen , 2009). This 
can, in turn, be useful to estimate a few system parameters 
(e.g., neutron star magnetic field; Ibragimov and Poutanen 
(2009)), and hence to make some of the methods described 
throughout this review more reliable. 

Finally, we note that other methods have been pro- 
posed in the literature to constrain the neutron star M — R 
space for AMPs (Bhattacharyya , 2001). Gilfanov et al. 
(1998) reported that the X-ray energy spectrum of SAX 
J1808. 4-3658 had remained stable, when the X-ray lumi- 
nosity varied by two orders of magnitude. Bhattacharyya 
(2001) (see also Burderi and King (1998); Li et al. (1999)) 
found the upper limits to the mass and the radius of the 
neutron star using this observational result, and with the 
following assumptions: (1) the neutron star magnetic field 
is dipolar; (2) there is no "intrinsic" pulse mechanism; and 
(3) the accretion flow is not centrifugally inhibited. 

3.3. Kilohertz Quasi-periodic Oscillations 

3.3.1. What are Kilohertz Quasi-periodic Oscillations? 

Neutron star LMXBs exhibit a zoo of timing features, 
both broad and narrow, in the Fourier power spectra (see 
van der Klis (2006) for an excellent review). Many of 
these features are correlated with each other and with the 
spectral states of the source (van der Klis (2006) and ref- 
erences therein). Among these features, high frequency 
quasi-periodic oscillations, or kilohertz (kHz) QPOs, are 
the most important ones to constrain the neutron star pa- 
rameters (see § 3.3.2). These QPOs often occur in a pair 
(Figure 11), and the twin peaks usually move together in 
the frequency range ~ 200 — 1200 Hz in correlation with 
the source state. The higher frequency QPO is known 
as the upper kHz QPO (frequency = v u ), and the lower 
frequency QPO is called the lower kHz QPO (frequency 
= v{). The kHz QPOs were first reported in van der Klis 
et al. (1996) and Strohmayer et al. (1996) (see also 
van der Klis (1998) for a historical account). As of now 
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these QPOs have been observed from about 30 neutron 
star LMXBs (see Altamirano et al. (2009a); van der Klis 
(2006)). 
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Figure 11: Power spectrum of a neutron star LMXB showing a pair of 
kHz QPOs (see § 3.3.1). The vertical dotted lines show the centroid 
frequencies of these QPOs. 

Although frequency (v), amplitude, quality factor (yj&v) 
and other properties of kHz QPOs vary from source to 
source, and even for a given source depending on its inten- 
sity and spectral state, several general trends of this timing 
feature could still be established. Here we mention some 
of them which may be useful to understand these QPOs. 
(1) For a given source, v correlates well with the X-ray 
luminosity L x on time scales of hours. But on longer time 
scales, and across sources, the v ~ L x correlation does not 
exist, and similar v values are observed over two orders of 
magnitude in L x (e.g., van der Klis (1997)). (2) It was 
thought that the frequency separation v u — vi (= A^) clus- 
ters around the neutron star spin frequency and the half 
of it (Wijnands ct al. (2003); van der Klis (2006) and 
references therein); although the recent works of Mendez 
and Bclloni (2007) and others have suggested that Av is 
consistent with no dependence on spin. (3) KHz QPO am- 
plitudes increase with photon energy, and in similar energy 
bands, the QPOs are weaker in the more luminous sources 
(Jonkcr ct al. (2001); van der Klis (2006) and references 
therein). (4) The quality factor of lower kHz QPOs in- 
creases with the QPO frequency at lower frequency values, 
and decreases at higher frequency values. But the quality 
factor of upper kHz QPOs increases with the frequency all 
the way to the highest detectable frequencies (Barret et al. 
, 2006). (5) KHz QPO frequencies are correlated with the 
frequencies (or characteristic frequencies) of many narrow 
and broad timing features from neutron star LMXBs (see 
Figure 2.9 of van der Klis (2006)). 

What causes the kHz QPOs is not known yet, although 
many models are available in the literature. A discussion 
of all these models is out of the scope of this review. Here, 
we will mention only a few models which will be used in 
§ 3.3.2 to demonstrate the plausible ways to measure neu- 
tron star parameters using kHz QPOs. The most promi- 



nent property of a kHz QPO is its high frequency. Most of 
the models have primarily attempted to explain this prop- 
erty. This frequency is usually believed to be caused by the 
frequencies related to the fast motions close to a neutron 
star, where the effects of general relativity is important. 
Let us first identify some of these frequencies, and write 
down their expressions for circular orbits in the equatorial 
plane for Kerr spacetime. 

(1) Orbital frequency around the neutron star 

vt = Ml+3(r 9 /rf/ 2 )-\ (9) 

where vk = \/GM/r 3 /2ir. Here, r is the radial distance 
from the center of the neutron star, r g = GM/c 2 , and the 
angular momentum parameter j = Jc/GM 2 , where J is 
the stellar angular momentum. 

(2) Radial epicyclic frequency (for infinitesimally eccentric 
orbits) 

Vr = ^(i - 6 (r fl /r) + 8j(r g /rf/ 2 - 3j 2 (r fl /r) 2 ) 1/2 - (10) 

(3) Vertical epicyclic frequency (for infinitesimally tilted 
orbits) 

u e = ^(1 - 4j(r s /r) 3 / 2 + Sf^/r) 2 ) 1 / 2 . (11) 

(4) Periastron precession frequency v pci i — — v r . 

(5) Nodal precession frequency f no dai = — vg. 

Apart from these, the neutron star spin frequency f sp in 
may also contribute to kHz QPOs. The radial profiles of 
the general relativistic frequencies are shown in the Fig- 
ure 12. This figure helps to understand any model based 
on these frequencies. It also clearly shows that the ra- 
dial epicyclic frequency cannot be v u or v\, because some 
kHz QPOs are observed with frequencies much higher than 
the maximum possible value of v r . In many models, kHz 
QPO frequencies are explained in terms of the above men- 
tioned general relativistic frequencies and f sp i n (or beating 
between any two of them) at some preferred radii. A pre- 
ferred radius may be determined by the disk flow struc- 
ture (depends on, for example, the accretion rate), the 
radius of the general relativistic innermost stable circular 
orbit (ISCO), or the commensurabilities among the fre- 
quencies which may cause resonance. For example, Miller 
et al. (1998) identified v u and v\ with and a beating fre- 
quency (fbeat) respectively at the "sonic" radius (r son ; c ). 
Lamb and Miller (2003) modified this model maintain- 
ing v u — fcj,(r aon i c ), and explaining v\ by a beat-frequency 
interaction between orbital motion at r son j c and an az- 
imuthal structure at the spin-resonance radius, defined by 
Z4 p in — V4, = vq. Stella and Vietri (1998, 1999) identi- 
fied v u with i/^ at the disk inner edge, and related v\ and 
vu (frequency of a low-frequency QPO; see van der Klis 
(2006)) with f per ; and ^ no dai respectively. Various authors 
(e.g., Kluzniak and Abramowicz (2001); Abramowicz and 
Kluzniak (2001)) made use of the fact that v r and vg 
at particular radii have simple integer ratios or other com- 
mensurabilities with each other or with v spin . They sug- 
gested that at these radii resonances may occur which can 
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show up as kHz QPOs. Besides, many other authors have 
used the general relativistic frequencies in various ways to 
explain the properties of kHz QPOs (e.g., Wijnands et al. 
(2003); Lee et al. (2004); Zhang (2004); Mukhopadhyay 
(2009)). 
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volve plasma motion in the strong gravitational field around 
the neutron star, and with one exception (photon bubbles; 
Klein et al. (1996a, b)) suggest that the kHz QPOs orig- 
inate in the disk (van der Klis , 2006). Moreover, most 
models identify one of the kHz QPO frequencies (usually 
v u ; but can also be v\) with the orbital motion at a pre- 
ferred disk radius (van der Klis (2006); see also § 3.3.1). 
If this is true, two reasonable conditions can constrain the 
M - R space (Miller et al. , 1998): 



R<r, 



(12) 



where r is the radius of the orbit associated with v u or v\ 
via equation 9; and 



nsco < r, 



(13) 



where nsco is the radius of the ISCO. This is because the 
first condition gives a mass-dependent upper limit on R 
via equation 9; and the second condition gives an upper 
limit on M: M < c 3 /{2n6 3/2 G^\ r ) (for Schwarzschild 
spacetime) . These constraints on M — R space are shown 
in Figure 13. 



Figure 12: This figure shows the radial profiles of various frequencies 
(color coded) of equatorial circular orbits in Kerr spacetime (see 
§ 3.3.1). Two angular momentum parameters (j = 0.0 (solid) and 
j = 0.3 (dotted)), and neutron star mass = 1.35M0 arc used. These 
frequencies may be useful to understand the kHz QPOs (see § 3.3.1). 

Apart from the observed frequencies of kHz QPOs, 
one also needs to explain the modulation and decoherence 
mechanisms. Some of the plausible methods have been 
discussed in van der Klis (2006) and references therein. 
Mendez (2006) proposed that, although the kHz QPO 
frequencies are plausibly determined by the characteristic 
disk frequencies (see the previous paragraph), the modu- 
lation mechanism is likely associated to the high energy 
spectral component (e.g., accretion disk corona, boundary 
layer between the disk and the neutron star, etc.). This is 
because the disk alone cannot explain the large observed 
amplitudes, especially at hard X-rays where the contribu- 
tion of the disk is small. Detection and measurement kHz 
QPOs above 20 — 30 keV may be able to resolve this issue. 
In addition, observations of sidebands, overtones, and very 
high frequency QPOs may be useful to identify the correct 
kHz QPO model (van der Klis , 2006; Bhattacharyya , 
2009). We will not discuss these aspects further. Rather 
we note that the fluid dynamical simulation corresponding 
to any successful kHz QPO model must naturally give rise 
to the observed frequencies and other properties. 

3.3.2. Kilohertz Quasi-periodic Oscillation Method 

Since no kHz QPO model can yet explain all the major 
aspects of this timing feature, currently it is not possi- 
ble to constrain the neutron star parameters using this 
QPOs with certainty. However, all proposed models in- 
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Figure 13: This figure shows the M — R space of neutron stars with 
the curves corresponding to a few representative EoS models (same 
as Figure 1). The green patch shows the allowed M — R space using 
equations 12 and 13, and an upper kHz QPO frequency of 1200 Hz 
(Miller et al. (1998); see § 3.3.2). This example figure demonstrates 
the potential of kHz QPOs to constrain the neutron star mass and 
radius, and hence the EoS models. 

The general relativistic frequencies depend on the neu- 
tron star parameters M and j, and many models use these 
frequencies to explain kHz QPOs (§ 3.3.1; see also van 
der Klis (2006) and references therein). Therefore, in 
a more general sense the identification of any of the kHz 
QPO frequencies with the beating, resonance, or any other 
combination of f^, v r , vg, is pci i, v no ded and i^pm nas the 
potential to constrain the neutron star parameter space. 
Besides, the identification ofvh with f no dai (§ 3.3.1) can be 
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useful to constrain the stellar moment of inertia /, because 
f nodal = ^ 2 v 2 Iv spin /Mc 2 (Stella and Vietri , 1998). Fi- 
nally, simultaneous observations of kHz QPOs and broad 
relativistic iron lines have the potential to measure neu- 
tron star parameters (see § 3.4.2). 

3-4- Broad Relativistic Iron Lines 

3.4-1- What is a Broad Relativistic Iron Line? 

A broad iron Ka spectral emission line near 6 keV is ob- 
served from many accreting supermassive and stellar-mass 
black hole systems (see Fabian et al. (2000); Reynolds and 
Nowak (2003); Miller (2007); and references therein). Al- 
though various models were proposed to explain this fea- 
ture (e.g., Misra and Kembhavi (1998); Titarchuk et al. 
(2003); Reeves et al (2004); Titarchuk et al. (2009)), 
the large width, the lower-energy wing and the asymme- 
try of such a line, as well as the observed associated "disk 
reflection" spectrum strongly suggest that this line origi- 
nates from the inner part of the accretion disk, where it 
is shaped by Doppler and relativistic effects (e.g., Tanaka 
et al. (1995); Fabian et al. (2000); Miller et al. (2001); 
Miniutti et al. (2007); Miller et al (2009)). This relativis- 
tic line are believed to be produced by the reflection of hard 
X-rays from the accretion disk. The hard X-ray source 
may be anything from an accretion disk corona (Fabian et 
al. , 2000) to the base of a jet (Markoff and Nowak , 2004; 
Markoff et al. , 2005). A given incident X-ray photon may 
be Compton scattered by free or bound electrons, or sub- 
ject to photoelectric absorption followed by either Auger 
de-excitation or fluorescent line emission (Fabian et al. , 
2000) . The strongest among the fluorescent spectral emis- 
sion lines is the one for the n = 2 — > n = 1 transition 
of the iron atom (or ion). This iron Ka line originates 
at an energy between 6.4 keV and 6.97 keV, depending 
on the iron ionization state (Reynolds and Nowak , 2003). 
The intrinsically narrow iron Ka line is shaped by the fol- 
lowing effects: (1) the Doppler effect due to the motion 
of the line-emitting matter in the inner part of the accre- 
tion disk broadens the line and makes it double-peaked; 
(2) The general relativistic light-bending effect also con- 
tributes to the formation of the two peaks; (3) the trans- 
verse Doppler shift and the gravitational redshift shift the 
line towards lower energies; and (4) the special relativis- 
tic beaming reduces the lower-energy peak and enhances 
the higher-energy peak making the line asymmetric (we 
encourage the readers to see the Figure 3 of Fabian et al. 
(2000)). Therefore, the relativistic iron line provides an 
ideal tool to probe the spacetime and the flow and struc- 
ture of accreted matter in the strong gravity region. The 
modelling of such a line can also be useful to constrain the 
radius of the ISCO (in the unit of black hole mass) , which 
can be used to measure the black hole spin parameter or 
angular momentum parameter (j = Jc/GM 2 = a/M; J 
and M are the black hole angular momentum and mass re- 
spectively; see, for example, Brenneman Reynolds (2006); 
Miller et al (2009)). This is because a/M has a one-to- 
one correspondence with the ratio of the ISCO radius to 



black hole mass for Kerr spacetime (see Figure la of Miller 
(2007)). 
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Figure 14: Broad relativistic iron Ka spectral line from the inner 
accretion disk of a neutron star LMXB Serpens X-l (Bhattacharyya 
and Strohmaycr , 2007b) . The x-axis shows the X-ray energy, and the 
y-axis gives the observed intensity in excess to the best-fit continuum 
spectral model. The data points clearly show a broad asymmetric 
spectral emission line, while the dotted profile is a relativistic model 
which fits the observed line well (§ 3.4). 

Similar to black hole systems, a broad iron line from 
neutron star LMXBs has also been known for many years 
(see, for example, Asai et al. (2000)). However, the in- 
ner accretion disk origin of this line could not be estab- 
lished until recently. This was because the characteristic 
asymmetry of the relativistic nature could not be detected 
due to the modest signal-to-noise ratio. In 2007, Bhat- 
tacharyya and Strohmayer (2007b), for the first time, es- 
tablished the inner accretion disk origin of the broad iron 
line from a neutron star LMXB (see Figure 14). This was 
done by analyzing the XMM-Newton data from Serpens 
X-l. Soon, Cackett et al. (2008a) confirmed this finding 
using the independent Suzaku data from the same source. 
These authors also reported the detection of relativistic 
iron lines from two other neutron star LMXBs: 4U 1820- 
30 and GX 349+2. As of now, the inner disk origin of 
broad iron line has been confirmed for ten neutron star 
LMXBs (Bhattacharyya and Strohmayer , 2007b; Cackett 
et al. , 2008a; Pandel et al. , 2008; D'Ai et al. , 2009; 
Cackett et al. , 2009a; Papitto et al. , 2009; Shaposhnikov 
et al. , 2009; Reis et al. , 2009; di Salvo et al. , 2009; Iaria 
et al. , 2009; Cackett et al. , 2009b). This opens up a new 
way to measure neutron star parameters, as described in 
§ 3.4.2. 

3-4-2. Broad Relativistic Iron Line Method 

The broad relativistic iron line can be used to constrain 
various neutron star parameters. This line originates from 
an annular region of the inner accretion disk (see § 3.4.1). 
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The ratio of the inner edge radius of this annulus (rj n!ann ) 
to the neutron star mass (M) can be constrained from the 
fitting of the iron line with a relativistic model (see, for 
example, Bhattacharyya and Strohmayer (2007b)). Since 
the disk inner edge radius has to be obviously greater 
than or equal to the neutron star radius, the measured 
nn,annC 2 /GM gives a hard upper limit on the neutron star 
Rc 2 /GM. This can be useful to constrain the stellar M—R 
space (see, for example, Figure 6). 

The relativistic iron line may also be used to constrain 
the neutron star angular momentum parameter, if ri D , ann c 2 /' 

< 6. Here is how. The accretion disk can extend up to 
the radius (nsco) of the ISCO, or the neutron star radius 
(R), whichever is greater (see Figure 1 of Bhattacharyya 
et al. (2000)). However, the disk may be truncated at a 
larger radius by the magnetic field and/or radiative pres- 
sure, and/or by any other means. Therefore, if r; n is the 
disk inner edge radius, then ri n . ann c 2 /GM > r- ln c 2 jGM > 
Rc 2 /GM and r injann c 2 /GM > r in c 2 /GM > n sco c 2 /GM. 
For a non-spinning neutron star, nscoc 2 /GM — 6, and for 
a spinning neutron star with a corotating disk, r\a,coc 2 / GM 

< 6. Therefore, if the best fit value of nn,ann c 2 /GM is less 
than 6, then that must be because of the effect of the stel- 
lar spin on the spacetime, and hence can, in principle, be 
useful to constrain the neutron star angular momentum 
parameter (see Figure la of Miller (2007)). Here we note 
that the Kerr spacetime has been assumed in the Figure 
la of Miller (2007), while the exterior spacetime of a spin- 
ning neutron star deviates from the Kerr geometry (Miller 
and Lamb , 1996). 

In § 3.3, we have mentioned that kHz QPOs are likely 
associated with the inner accretion disk. If this is true, 
then the relativistic iron line and the kHz QPOs are possi- 
bly physically related, and their simultaneous observation 
and joint analysis may, not only shed light on the physics 
of these features, but also be useful to constrain the neu- 
tron star parameters. For example, Miller et al. (1998) 
suggested that the upper kHz QPO frequency (v n ) is close 
to the orbital frequency at the inner edge of the optically 
thick emission, and thus from a radius close to the r; njann 
inferred from the iron line. If this is true, then the neu- 
tron star mass (M) can be measured from the relation 

u u = (l/27r)y f GM /r? n ann (assuming Schwarzschild space- 
time; Cackett et al. (2008a)). 

Broad relativistic iron line method has emerged as a 
way to constrain the neutron star parameters only re- 
cently. Although this method has a great potential, it 
suffers from systematics similar to other methods. These 
systematic uncertainties may appear due to (1) the lack 
of knowledge of various parameters, such as the observer's 
inclination angle, the nature of the hard X-ray emission 
incident upon the disk, etc; and (2) the inability to model 
the continuum spectrum accurately. Here we will not dis- 
cuss these uncertainties in detail. Rather, we will only 
note that the incident hard X-ray emission in the neutron 
star LMXBs is currently believed to be originated from 



the boundary layer between the inner accretion disk and 
the neutron star surface (Cackett et al. , 2009b), although 
more studies are required to verify this. 

3.5. Quiescent Emissions 

3.5.1. What is Quiescent Emission? 

Many neutron star LMXBs have two distinctly differ- 
ent intensity states: (1) the quiescent state in which the 
source remains for months to years; and (2) the outburst 
state which continues for weeks to months (sometimes for 
fears; e.g., for the source KS 1731-260). These LMXBs 
are called transients. The luminosity of the source typi- 
cally increases by several orders of magnitude as it evolves 
from the quiescent state into an outburst. Figure 1 of 
Bhattacharyya et al. (2006c) shows the long term in- 
tensity profile of such a transient 1A 1744-361, exhibiting 
three outbursts. The transient nature of these LMXBs is 
normally attributed to an accretion disk instability, which 
causes high accretion rates for certain periods of time (out- 
bursts) , and almost no accretion during other times (quies- 
cence; see King (2001) for a review). Therefore during the 
quiescent state, the observed X-ray emission is expected 
to primarily originate from the neutron star surface. This 
emission can be detected and measured with high sensitiv- 
ity imaging instruments on board Chandra, XMM-Newton, 
etc. (e.g., Wijnands et al. (2002); Cackett et al. (2008b); 
Heinke et al. (2009)). 

3.5.2. Quiescent Emission Method 

As mentioned in § 3.5.1, the emission from a transient 
neutron star LMXB in quiescence primarily comes from 
the neutron star's atmosphere. This atmosphere should be 
composed of pure hydrogen, and devoid of any X-ray spec- 
tral line (see, for example, van Kerkwijk (2004)). More- 
over, the relatively low magnetic fields of neutron stars in 
LMXBs (unlike isolated neutron stars, and stars in dou- 
ble neutron star binaries and HMXBs) make the modeling 
simpler. Therefore, assuming a blackbody emission from 
the neutron star surface in a quiescent LMXB, one can 
constrain the stellar radius using equations 2 and 4. As 
mentioned in § 3.1.3, there may be effects of systematics in 
the estimated radius. During quiescence, the entire neu- 
tron star surface is expected to emit uniformly, because the 
stellar magnetic field is low, and there is no obvious effect 
which can create a significant asymmetry. Therefore, The 
systematic uncertainty #1 (§ 3.1.3) may not be present. 
The source distance, which appears in equation 2, may be 
known from a PRE burst from the source, or in case the 
source is in a globular cluster (see uncertainty #2; § 3.1.3). 
The surface gravitational redshift 1 + z, which appears in 
equation 4 and causes the uncertainty #3 (§ 3.1.3), may be 
determined from an independent measurement of the neu- 
tron star radius-to-mass ratio (see, for example, § 3.1.4; 
§ 3.1.6; § 3.4.2). But, even if this ratio is not known, 
the neutron star M — R space can still be constrained as 
shown in Figure 4. The color factor /, which appears in 
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equation 4 and causes the uncertainty #4 (§ 3.1.3), may 
be determined by theoretical calculations. However, if a 
theoretical atmospheric spectrum for hydrogen is used for 
fitting (instead of a blackbody spectrum), then a sepa- 
rate computation of / is not required. In fact, hydrogen 
atmospheric spectrum is usually used to model the quies- 
cent emission, which can constrain the inferred radius Roo 
of the neutron star well (equations 2 & 4; Gendre et al. 

(2003) ; Webb and Barret (2007); Guillot et al. (2009)). 
Apart from the above systematics, three observational 

aspects make the quiescent emission method less reliable. 
(1) The interstellar absorption cuts off a large part of 
the energy spectrum for some sources (van Kerkwijk , 
2004). This makes the correct spectral fitting difficult. 
One needs to measure the interstellar absorption indepen- 
dently (for example, by simultaneous UV observation) to 
solve this problem. (2) The thermal component of the 
energy spectrum can vary in a timescale of months (Rut- 
ledge et al. , 2002). Moreover, in addition to the thermal 
emission, sometimes a non-thermal component is found 
(see, for example, Rutledge et al. (2002)). This non- 
thermal emission may be due to a residual accretion (van 
Kerkwijk , 2004). Detailed observations of many quies- 
cent LMXBs with XMM-Newton, Chandra and future X- 
ray instruments will be required to interpret the quiescent 
emission X-ray spectrum correctly, and hence to increase 
the reliability of the method. (3) Episodic higher accre- 
tion rates have been inferred during the quiescent emission 
(Kuulkers et al. , 2009). However, such high accretion is 
usually followed by a thermonuclear burst (Kuulkers et al. 
, 2009). Therefore if such a burst is detected, then the data 
immediately prior to it can be excluded from the neutron 
star parameter measurement analysis. 

Finally we note that a neutron star heats up during 
accretion, and cools down when the system goes into the 
quiescence. The cooling curve of the neutron star can be 
useful to understand its properties (e.g., Wijnands et al. 

(2004) ). Since repeated outbursts can be observed from a 
given transient system, the neutron star can heat up and 
cool down repeatedly, providing a useful tool to study the 
stellar core and crust. 

3.6. Mass Measurement: Binary Orbital Motion Method 

The mass M of a neutron star affects the motion of 
the companion star via gravitation, and hence this mass 
may be constrained from the optical observations of the 
signatures of binary orbital motions. For example, the 
radial velocity (RV; i.e., the orbital velocity component 
along the line of sight) curve of the companion, based on, 
say, its absorption line spectra, can be useful to measure 
the binary orbital period (P or b) and the neutron star mass 
function: 

/(M) _ _ Mj^ (14) 



is the companion mass. Note that, since q > and 



(l + <?) 2 ' 



sini < 1, f(M) gives the lower limit of the neutron star 
mass. Furthermore, with the estimated i and q, M can be 
constrained from equation 14. The inclination angle i can 
be estimated from the observed light curves, and/or from 
X-ray dips/eclipses that may be detected from the LMXB 
(e.g., White and Swank (1982); Parmar et al. (1986)). 
For example, observations of both dips and total eclipses 
may imply i w 75° - 80° (Frank et al. , 1987). Dips and 
eclipses also provide the binary orbital period (e.g., White 
and Swank (1982); Parmar et al. (1986)). The mass ratio 
q of equation 14 can be estimated by measuring the rota- 
tional (i.e., spin- induced) broadening of the spectral lines 
originated from the companion star. This is because the 
spin-related linear speed (o comp ) of the companion surface 
can be expressed in terms of q using the following formula 
(Wade and Home , 1988): 



,sini = 0.462X comp(? 1 / 3 (l + qf' z . 



(15) 



However, if an estimate of q is not available, then one needs 
to measure the companion star mass function: 



f(M c , 



, _ -Porb ^ NS _ M comp sin 3 j 

2nG _ (l + [l/g])2 



(16) 



Here K comp is the amplitude of the companion star RV, 
i is the observer's inclination angle, q = M comp /M, and 



in order to constrain the neutron star mass. Here .Kns is 
the amplitude of the RV of the neutron-star. 

Unfortunately, the optical emission of a neutron star 
LMXB is heavily dominated by the accretion disk radi- 
ation for most sources. This optical radiation originates 
from the reprocessing of the X-rays in the outer accre- 
tion disk. As a result, the above mentioned binary orbital 
motion method cannot be effectively used for the neutron 
star LMXBs, except for a few sources for which the com- 
panion star can be detected, and its relevant properties 
can be measured sufficiently well (e.g., Orosz and Kuulk- 
ers (1999); Jonker et al. (2005)). An example of such 
a system is Cyg X-2, because s=s 70% of its optical emis- 
sion originates from the companion star (Elcbcrt et al. , 
2009). The transient neutron star LMXBs in their qui- 
escent state might therefore be ideal to apply the above 
mentioned method (see § 3.5). But, at a distance of a few 
kpc the low-mass companion stars are usually too faint to 
observe, again making the binary orbital motion method 
ineffective. 

In 2002, Steeghs and Casares (2002) discovered nar- 
row high-excitation optical emission lines from the neu- 
tron star LMXB Sco X-l. The strongest lines were in the 
Bowen range (4630A— 4660A), that primarily consisted of 
a blend of N III and C III lines (see Figure 15). The 
N III lines originate from a UV fluorescence process, and 
the C III lines are due to photo-ionization and subsequent 
recombination (McClintock et al. , 1975). These lines 
are believed to originate from the irradiated surface of 
the companion star, mainly because of the following rea- 
sons: (1) the Doppler tomography of Bowen blend emis- 
sion lines revealed a bright spot of emission (see Fig. 4 of 
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Figure 15: Average normalised optical spectrum of Sco X-l in the 
rest frame of the companion star. The Bowen blend N III and C 
III emission lines from the irradiated surface of the companion star, 
and the He II emission line are shown (§ 3.6; figure courtesy: Danny 
Steeghs; Steeghs and Casares (2002)). 

Steeghs and Casares (2002)); and (2) phase-resolved spec- 
troscopy showed that these narrow lines moved roughly 
in anti-phase with the neutron-star/accretion-disk com- 
ponent (see Figure 16). As a result these lines led to 
the first estimate of K comp for a neutron star LMXB for 
which the optical emission primarily originates from the 
accretion disk. Therefore, the discovery of Steeghs and 
Casares (2002) and the subsequent applications of the 
above mentioned binary orbital motion method using the 
Bowen blend emission lines for several sources have es- 
tablished a new technique to constrain the neutron star 
masses (see Cornelisse et al. (2008) and references therein; 
see also Muhoz-Darias et al. (2009)). 

This new technique, although is very promising, suffers 
from systematics mainly because of the unknown values of 
various source parameters. For example, i and q of equa- 
tion 14 and -Kns of equation 16 are not known for many 
sources. Besides, K comp measured from the narrow Bowen 
blend emission lines is only a lower-limit to the true K comp , 
because these lines arise on the irradiated side of the com- 
panion star which does not correspond to the center of 
mass of the companion. This shift of the measured -Kcomp 
from the true K comp can be accounted for by the "K- 
correction", which depends on i, q and the opening-angle 
of the accretion disk (Muiioz-Darias et al. (2005); see also 
de Jong et al. (1996)). Among these source parameters, 
the uncertainty on i and the opening-angle of the accretion 
disk could be reduced by modeling the accretion disk, and 
q could be obtained by clearly resolving the narrow lines 
to determine the rotational broadening. Another system- 
atic uncertainty comes from the poorly measured systemic 
velocity, which is an important parameter to do Doppler 
tomography (e.g., Cornelisse et al. (2009)). Higher resolu- 
tion data can be useful to reduce this uncertainty. Finally, 
the Bowen blend emission lines may also primarily orig- 
inate from the gas-stream/accretion-disk impact region, 
rather than from the irradiated surface of the companion 
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Figure 16: Upper panel: Radial velocity (RV; in km/s) curve of 
the narrow Bowen blend emission lines (see Figure 15) from Sco X- 
1. Data are repeated over two cycles, and the best sinusoidal fit is 
plotted as a dash-dot line. Lower panel: Data and sinusoidal fit of 
the RV curve of the He II emission line (see Figure 15) from Sco X-l. 
Since the Bowen blend emission lines originate from the irradiated 
surface of the companion star, and the He II line comes from near 
the neutron star, the RV curve of the latter is naturally roughly in 
anti-phase with that of the former (§ 3.6; figure courtesy: Danny 
Steeghs; Steeghs and Casares (2002)). 

star, for some neutron star LMXBs (Elcbcrt et al. , 2009). 

4. Why are Low-mass X-ray Binaries Useful? 

As mentioned in § 1 (see also Figure 1), one needs to 
measure at least three independent parameters of the same 
neutron star in order to constrain the EoS models of the 
supranuclear core matter. Because of a number of system- 
atics (mentioned throughout this review), it is possible 
only if several measurement methods can be used for a 
given neutron star. This may be achieved for neutron star 
LMXB systems, because many of them exhibit several of 
the following features: thermonuclear bursts without oscil- 
lations (suitable for continuum spectrum method), burst 
oscillations, PRE bursts, mHz QPOs, accretion-powered 
millisecond period pulsations, kHz QPOs, broad relativis- 
tic iron lines, quiescent emission, Bowen blend emission 
lines, etc. For example, SAX J1808. 4-3658 has shown all 
these features except mHz QPOs, and 4U 1636-536 has ex- 
hibited all these features except accretion-powered pulsa- 
tions and quiescent emission. Therefore, Bayesian method 
may be used to strongly constrain the neutron star param- 
eter values, which do not change from one observation to 
another (see § 3.1.8 for a discussion). This is why LMXBs 
are useful systems to measure the neutron star parame- 
ters. In Table 1 we give examples of a few sources for 
which various methods mentioned in this review have been 
attempted, and in Figure 17 we demonstrate how two of 
these methods can be jointly used to tightly constrain the 
M — R space. 



21 



5. Summary and Future Prospects 
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Figure 17: This figure shows the M — R space of neutron stars with 
the curves corresponding to a few representative EoS models (same 
as Figure 1). Miller et al. (1999) mentioned that the highest kHz 
QPO frequency observed from the neutron star LMXB 4U 1636-536 
was 1220 Hz, which gives an allowed M — R region for this source 
according to the equations 12 and 13 (compare this figure with Fig- 
ure 13). In addition, Nath et al. (2002) found GM/Rc 2 < 0.163 
for 4U 1636-536 by modeling the burst oscillations. The green patch 
shows the intesection of these two allowed regions, which should be 
the allowed M — R space for 4U 1636-536 using the kilohertz quasi- 
periodic oscillation method (§ 3.3.2) and the burst oscillation method 
(§ 3.1.6). Therefore, this figure demonstrates how various methods 
described throughout this review can be jointly used to tightly con- 
strain the M — R space. However, note that the particular allowed 
M — R space shown here may not be very reliable, because the ac- 
tual origin of kHz QPOs is not yet known, and the burst oscillation 
measurement of M/R used here suffered from systematics. 

There is another reason that makes neutron star LMXBs 
useful. In order to understand the nature of the supranu- 
clear core matter extremely well, we will need to observa- 
tionally trace out the stellar M — R curve (see, for exam- 
ple, Ozel and Psaltis (2009); see also Figure 1). Although, 
currently we cannot do it, with the technological advance- 
ment it may be doable in the future. But this will be pos- 
sible only if neutron stars have a sufficiently large range 
of mass values. Since the neutron stars in LMXBs can 
accrete mass plausibly up to 0.7M Q (e.g., van den Heuvel 
and Bitzaraki (1995)), they should occupy a larger range 
in the mass space than the isolated neutron stars, and the 
stars in double neutron star binaries. These latter neutron 
stars are formed in HMXBs, and hence a limited amount 
of matter can be accreted because of the short life of the 
massive companion. In fact, Thorsett and Chakrabarty 
(1999) found neutron star mass M = 1.35 ± 0.04M© for a 
sample of double neutron star binaries. Besides, LMXBs 
may contain neutron stars massive enough to rule out the 
softer EoS models (see § 1). Therefore, the LMXB sys- 
tems may be intrinsically more useful than the isolated 
stars, and stars in double neutron star binary systems. 



In this review, we have described a few methods to con- 
strain the neutron star parameters. We have chosen the 
methods which utilize various properties of thermonuclear 
bursts, accretion-powered millisecond pulsations, kHz QPOs, 
broad relativistic iron lines, quiescent emissions, and bi- 
nary orbital motions. Some of these phenomena provide 
more than one method. For example, thermonuclear burst 
alone provides five methods. Throughout the review we 
have stressed that various systematic uncertainties hinder 
the accurate measurement of the neutron star parameters, 
and a joint analysis of several phenomena has a great po- 
tential to reduce these uncertainties. The large amount of 
archival neutron star LMXB data accumulated with the 
current and past X-ray space missions can be used for 
such a joint analysis, and this analysis can be further im- 
proved with the advanced statistical and other techniques 
(e.g., Noble and Nowak (2008)). These techniques may in- 
clude pulsed phase spectroscopy, higher order timing anal- 
ysis, and a higher capability spectral fitting recipe, which 
can, for example, handle the fitting of burst continuum 
spectra with a large number of realistic atmospheric mod- 
els. However, such techniques should be complemented 
with theoretical investigations of the various phenomena. 
For example, astrophysicists and nuclear physicists should 
work more cohesively to understand the various properties 
of thermonuclear bursts. Numerical study of the inner ac- 
cretion disk in order to jointly model the kHz QPOs and 
the broad relativistic iron lines will also be useful. Besides, 
processes of the channelled flow and the radiative transfer 
in AMPs should also be numerically investigated. 

In addition to the archival data analysis complemented 
with advanced techniques and theoretical studies, the pro- 
posed X-ray space missions also hold a great promise. The 
Indian multiwavelength astronomy space mission Astrosat 
is planned to be launched in 2010, and will observe simul- 
taneously in a wide energy range (optical to hard X-ray 
of 100 keV). Apart from this unprecedented capability, its 
LAXPC instrument should be able to detect and measure 
the high frequency timing features, such as kHz QPOs and 
accretion-powered pulsations, for the first time in hard X- 
rays (say, up to w 50 keV). This will be useful, among other 
things, to discriminate among various kHz QPO modu- 
lation mechanism models. The Japanese X-ray satellite 
Astro-H is scheduled to be launched in 2013. This mis- 
sion will have an unprecedented imaging capability in the 
0.3 — 80 keV range, and will increase the reliability of the 
neutron star parameter measurement methods. The In- 
ternational X-ray Observatory (7XO), jointly proposed by 
NASA, ESA and JAXA, will have several highly sensitive 
instruments which will take the X-ray study of neutron 
star LMXBs to a higher level. The launch of this space 
mission is planned for 2021. The proposed "High Timing 
Resolution Spectrometer" (HTRS) of 1XO will be able to 
observe sources with fluxes of 10 6 counts per second in 
the 0.3 — 10 keV band without performance degradation, 
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Figure 18: The residuals (x 2 deviations) for a simulated iron Lya 
absorption line from the neutron star surface during a thermonuclear 
X-ray burst. The assumed values are the following: surface gravita- 
tional rcdshift 1 + z = 1.35, line equivalent width = 12 eV, exposure 
= 100 s, burst blackbody temperature = 1.7 keV and X-ray flux = 
4 X 10 — 8 ergs s -1 cm -2 . The simulation was done for the "High 
Timing Resolution Spectrometer" (HTRS; previously proposed for 
the X-Ray Evolving Universe Spectrometer (XEUS), but now pro- 
posed for the International X-ray Observatory (IXO)). This figure 
shows that even a short observation of thermonuclear X-ray bursts 
with HTRS may lead to a significant detection of surface atomic 
spectral lines (§ 3.1.4; figure courtesy: Didier Barret; Barret et al. 
(2008)). 

while providing good spectral resolution and excellent time 
resolution (e.g., Barret et al. (2008)). This instrument 
should, therefore, be able to reduce the systematics sig- 
nificantly, and measure the parameters of several neutron 
stars accurately. In Figures 18, 19 and 20, we show the 
simulated capabilities of HTRS. The Advanced X-ray Tim- 
ing Array (AXTAR), which is an X-ray observatory mis- 
sion concept currently under study, will provide a better 
sensitivity than IXO for timing observations of accreting 
neutron stars (see Ray et al. (2009)). Finally, future 
large optical/IR telescopes, for example, the proposed Ex- 
tremely Large Telescope (Spyromilio et al. , 2008), Thirty 
Meter Telescope (Nelson and Sanders , 2008), can be very 
useful for the binary orbital motion method, especially by 
detecting the companion stars of transient neutron star 
LMXBs in the quiescent state. 
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Figure 19: The residuals (x 2 deviations) for simulated broad rel- 
ativistic iron lines at two different inner accretion disk radii. The 
simulation was done for 1000 s of HTRS exposure. This figure shows 
that a short HTRS observation of relativistic lines may significantly 
reveal their detailed structures (§ 3.4; figure courtesy: Didier Barret; 
Barret et al. (2008)). 
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Tabic 1: A few examples of neutron star LMXBs for which various parameter measurement methods have been attempted. Note that this 
list is not at all exhaustive, 



No. 


Method (Section)" 


Source (Reference) b 


1 


Burst Continuum Spectrum 
Method (§ 3.1.3) 


EXO 0748-676 (Ozcl , 2006); 

4U 1746-37 (Sztajno et al. , 1987); 

4U 1636-536 (Fujimoto and Taam , 1986) 


2 


Burst Spectral Line Method 
(§ 3.1.4) 


EXO 0748-676 (Cottam et al. , 2002; Ozcl , 2006) 


3 


Photospheric Radius Expansion 
Burst Method (§ 3.1.5) 


EXO 0748-676 (Ozel , 2006); 

4U 1636-536 (Fujimoto and Taam , 1986); 

X2127+119 (Smale , 1998); 

Cygnus X-2 (Titarchuk and Shaposhnikov , 2002) 


4 


Burst Oscillation Method (§ 3.1.6) 


4U 1636-536 (Nath et al. , 2002); 

XTE J1814-338 (Bhattacharyya et al. , 2005) 


5 


Millihcrtz Quasi-periodic 
Oscillation Method (§ 3.1.7) 


Not attempted yet. 


6 


Accretion-powered Millisecond 
Pulsation Method (§ 3.2.2) 


SAX J1808.4-3658 (Poutanen and Gicrlihski , 2003; Leahy et al. , 2008); 
XTE J1814-338 (Leahy et al. , 2009) 


7 


Kilohcrtz Quasi-periodic 
Oscillation Method (§ 3.3.2) 


4U 1636-536 (Miller et al. , 1999); 

Sco X-l (Abramowicz et al. , 2003; Stella and Vietri , 1999); 
GX 5-1 (Stella and Vietri , 1998); 
GX 17+2 (Stella and Vietri , 1998); 
SAX J1808.4-3658 (Zhang , 2009) 


8 


Broad Relativistic Iron Line 
Method (§ 3.4.2) 


Serpens X-l (Bhattacharyya and Strohmayer , 2007b); 
4U 1820-30 (Cackett et al. , 2008a); 
GX 349+2 (Cackett et al. , 2008a); 
Cygnus X-2 (Cackett et al. , 2009b) 


9 


Quiescent Emission Method 
(§ 3.5.2) 


Sources in the globular clusters: omega Cen, M13, NGC 2808, NGC 6304 
(Webb and Barret , 2007; Guillot et al. , 2009) 


10 


Binary Orbital Motion Method 
(§ 3.6) 


EXO 0748-676 (Munoz-Darias et al. , 2009); 
4U 1822-37 (Muhoz-Darias et al. , 2005); 
4U 1254-69 (Barnes et al. , 2007) 



"Neutron star parameter measurement methods and the sections of this review in which these methods have been described. 
'A few example neutron star LMXBs for which the methods have been attempted; and the corresponding example references. 
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